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Vju.  MisTitAtr  j 

s  Theoretica 1  and  experimental  studies  were  performed  to  obtain  fundamental  understanding 
of  initiation  and  propagation  of  low-velocity  detonation  (LVD)  in  confined  energetic 
liquids,  LVD  was  modeled  as  a  reactive  shock  propagation  in  a  liquid  cavitatcd  by 
interactions  with  the  container  wall.  Reaction  required  to  support.  Lhe  wave  in  this 
cavitation  model  is  initiated  by  bubble  collapse  in  the  compressed  cavitation  field. 

2 ■•:)  computer  studies  of  a  shock  propagating  in  a  contained  liquid  showed  precursor 
wall  waves  produce  regions  of  tension  sufficient  to  cavitate  the  liquid  before  it  is 
compressed  by  the  main  shock.  Simulation  of  hot-spot-initiated  reaction  by  addition 
of  f  nerrv  '  ehind  the  shock  substantiated  the  hypothesis  of  partial  energy  release  in 
l.Vl).  'Methods  of  calculating  LVD  parameters  were  formulated  with  the  bydrothcrmodynamit'S 
of  Ch?  cavitation  model,  and  a  complete  equation  of:  state  for  nitromethano  was 
introduced  into  the  TIGER  code  to  permit  detonation  calculations  for  incomplete 
reaction.  Calculations  showed  about  20  percent  of  liquid  is  required  Lo  support  LVD. 
Dynamics  of  collapsing  va(.or  bubbles  were  analyzed  to  study  hot-spot  reaction 
intiation  and  conditions  for  LVD  initiation.  Initial  rate  of  temperature  increase  in 
babbles  with  heat  ano  muss  transfer  was  found  to  bo  10  Lo  20  times  less  than  in 
bubbles  without  transfer  at  the  bubble  wall.  Pressure -time  histories  obtained  in 
confined  shocked  inonopi opellants  shewed  shock  decoy ,  initiation  of  LVD ,  and  Lire 
initiation  of  high-velocity  detonation .  LVD  peak  pres  .ores  were  4-  to  7~kbar,  and  LVD 
velocities  were  about  2.0  nm/piset  in  et.hyl  nitrate  and  3.0  rma/^isec  in  FEFO. 
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ABSTRACT 


Theoretical  and  experimental  studies  were  performed  to  obtain  a 
fundamental  understanding  of  the  initiation  and  propagation  of  low-velocity 
detonation  (LVD)  in  confined  energetic  liquids,  LVD  was  modeled  as  a 
reactive  shock  propagating  in  a  liquid  cavita  ,:.i  by  interactions  with  the 
container  wall.  The  reaction  required  to  str  ;rt  the  wave  in  this  cavi¬ 
tation  model  is  initiated  by  bubble  collapse  in  the  compressed  cavitation 
field. 

Two-dimensional  computer  studies  of  a  shock  propagating  in  a  liquid 
contained  in  a  cylindrical  steel  tube  showed  that  precursor  wall  waves 
produce  regions  of  tension  sufficient  to  cavitate  the  liquid  before  it  is 
compressed  by  the  main  shock.  Simulation  of  hot-spot-initiated  reaction 
by  addition  of  energy  behind  the  shock  substantiated  the  hypothesis  of 
partial  energy  release  in  LVD.  Methods  of  calculating  LVD  parameters 
were  formulated  with  the  hydrothermodynamics  of  the  cavitation  model, 
and  a  complete  equation  of  state  for  nitromethane  was  introduced  into  the 
TIGER  code  so  that  detonation  calculations  could  be  made  for  incomplete 
reaction.  The  calculations  showed  that  only  about  20  percent  of  the 
liquid  is  required  to  react  to  support  LVD.  The  dynamics  of  a  collapsing 
vapor  bubble  were  analyzed  to  study  the  hot-spot  initiation  of  reaction 
and  the  conditions  for  initiation  of  LVD.  The  initial  rate  of  temperature 
increase  in  a  bubble  with  heat  and  mass  transfer  was  shown  in  general  to 
be  10  to  20  times  less  than  that  in  a  bubble  without  transfer  at  the  wall. 

Experimental  studios  resulted  in  measurements  of  pressure  in  LVD. 
Pressure-time  histories  in  shocked  monopropellants  contained  in  cylindrical 
steel  tubes  were  recorded  by  piezoresistant  ytterbium  in-material  stress 
gages.  Gage  records  obtained  by  varying  the  monopropellant,  the  wall 
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thickness-to-diameter  ratio  of  the  tube,  and  the  initial  shock  pressure 
showed  shock  decay,  initiation  of  LVD,  and  the  initiation  of  high- 
velocity  detonation.  The  peak  pressures  recorded  by  LVD  were  in  the 
4-  to  7-,kbar  range,  and  the  LVD  velocities  of  propagation  were  about 
2.0  mm/psec  in  ethyl  nitrate  and  3.0  mm/psec  in  FEPO. 
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I  INTRODUCTION 


This  report  presents  the  results  of  a  three-year  theoretical  and 
experimental  program  whose  objective  was  to  obtain  a  fundamental  under¬ 
standing  of  low-velocity  detonation  (LVD)  phenomena  in  liquid  propellants, 
fuels,  and  explosives.  Low-velocity  detonation  occurs  in  confined,  highly 
energetic  liquids.  It  can  be  initiated  by  much  smaller  shock  impulses 
than  can  high-velocity  detonation  (HVD) ,  but  can  cause  considerable 
damage.  Present  hydrothermodynamic  detonation  theory  predicts  the  exis¬ 
tence  of  only  HVD  for  any  particular  thermodynamic  state  of  the  explosive. 
The  existence  of  LVD  therefore  represents  an  anomaly  heretofore  unexplained 
in  detonation  theory.  A  basic  understanding  of  LVD  is  of  great  practical 
significance  in  establishing  safety  requirements  for  the  handling  of 
liquid  fuels,  propellants,  and  oxplosives. 

The  following  tasks  were  accomplished  during  the  research  program: 

(1)  A  complete  literature  review  and  a  reevaluation  of  existing 
LVD  data  were  accomplished  to  identify  the  Individual  mecha¬ 
nisms  contributing  to  LVD  initiation  and  propagation.  The 
experimentally  observed  features  of  LVD  were  identified  and 
used  to  evaluate  the  various  models  of  LVD  that  have  been 
proposed.  A  model  that  offered  the  most  plausible  explanation 
of  the  phenomenon  and  that  was  consistent  with  the  experimental 
observation  was  adopted  as  a  starting  point  for  the  subsequent 
work.  This  model  was  called  the  cavitation  model  of  LVD. 

(2)  A  preliminary  theoretical  analysis  of  the  interaction  of  the 
confinement  system  with  the  liquid  explosive  was  carried  out. 

The  studies  established  that  precursor  wall  waves  capable  of 
causing  cavitation  in  the  liquid  ahead  of  the  main  liquid 
shock  are  possible. 

(3)  The  physical  and  chemical  properties  of  liquids  known  to  undergo 
LVD,  as  well  as  those  in  which  LVD  was  not  observed,  were 
examined.  The  results  of  this  task  suggested  the  significance 
of  cavitation  in  the  LVD  phenomenon. 
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(4)  Detailed  computer  calculations  utilizing  the  cavitation  model 
of  LVD  were  carried  out  to  study  the  cavitation  production 
mechanism,  the  effects  of  partial  energy  release,  the  initia¬ 
tion  of  LVD,  and  the  possible  transition  of  LVD  to  HVD, 

(5)  A  hydrothermodynamic  treatment  of  the  cavitation  model  of  LVD 
was  undertaken  in  an  attempt  to  develop  a  satisfactory  method 
for  calculating  LVD  propagation  velocity.  Interaction  of  the 
shock  and  the  .cavitation  field  was  accounted  for  simply  by 
ordering  the  time  scales  for  the  attainment  of  mechanical, 
thermal,  and  chemical  equilibrium.  Since  results  based  on 
the  attainment  of  chemical  equilibrium  did  not  agree  with 
experiment,  the  case  of  incomplete  reaction  was  considered. 

(6)  Characteristic  times  for  bubble  collapse  ■,  heat  transfer  within 
a  bubble,  and  heat  transfer  to  the  surrounding  liquid  were 
compared  to  assess  the  relative  rates  of  possible  processes 
occurring  in  the  compressed  cavitation  field.  The  initial 
collapse  of  a  vapor  bubble  was  treated  analytically  in  an 
attempt  to  obtain  a  more  detailed  account  of  bubble  dynamics 
and  to  determine  conditions  for  the  onset  of  reaction  and  the 
initiation  of  LVD. 

(7)  Experiments  using  metriol  trinitrate  (TMETN)  were  carried  out 
to  examine  the  effects  of  confinement  upon  LVD  initiation  and 
propagation. 

(8)  Experiments  were  carried  out  to  measure  the  stress  field  ahead 
of  and  during  the  passage  of  the  liquid  shock.  Detonation 
velocity  measurements  as  well  as  damage  criteria  were  used 

to  verify  the  observed  LVD  phenomenon. 

As  a  result  of  the  above  studies,  a  quantitative  model  of  LVD  has 
been  developed  that  is  consistent  with  the  experimentally  observed  proper¬ 
ties  of  LVD.  Further,  experimental  verification  of  some  of  the  more 
important  quantitative  predictions  of  the  model  has  been  carried  out. 

The  results  of  this  work,  intended  for  publication  in  technical 
journals,  are  presented  in  the  manuscripts  attached  as  Appendices  to 
this  report,  and  a  complete  list  of  references  of  the  work  performed 
during  the  present  research  program  is  given  in  the  Bibliography. 
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II  ASSESSMENT  OF  PAST  WORK 


For  most  practical  applications,  condensed  explosives  are  required 
to  behave  reproducibly ,  to  undergo  HVD,  and  to  produce  pressures  in  the 
100-kbar  region.  As  a  result,  investigations  of  detonation  phenomena 
have  been  focused  largely  on  HVD — rather  than  on  LVD,  which  produces 
pressures  in  the  10-kbar  region.  Consequently,  at  the  beginning  of  this 
research  program,  the  conditions  for  initiation  and  propagation  of  LVD 
were  not  as  well  established  as  were  those  for  initiation  and  propaga¬ 
tion  of  HVD.  This  disparity  in  our  understanding  was  accentuated  by 
the  fact  that  previous  experimental  studies  of  LVD  had  produced  only 
limited  results,  generally  of  a  qualitative  nature.  For  instance,  it  was 
clear  that  a  reactive  wave  traveling  at  about  sonic  speed  in  virgin 
explosive  was  associated  with  stable  LVD,  that  LVD  could.be  initiated 
more  easily  than  HVD,  and  that  initiation  depended  on  the  properties  and 
geometry  of  the  container  holding  the  explosive.  There  was,  however,  no 

physical  model  for  LVD  corresponding  to  the  Zeldovich-vonNeumann-Doering 

1-3* 

model  for  HVD.  Neither  were  there  criteria  for  calculating  detona¬ 
tion  velocity  from  the  properties  of  the  explosives  corresponding  to  the 
Chapman-Jouguet  (C-J)  hypothesis  nor  a  satisfactory  model  for  treating 
the  coupling  between  the  container  walls  and  the  explosive  in  the  LIT) 
process. 

During  the  initial  stages  of  this  program,  a  comprehensive  survey 
and  evaluation  was  made  of  existing  experimental  data  and  theoretical 
models  pertaining  to  LVD.  The  results  of  this  study  are  given  in  detail 


* 

References  are  listed  at  the  end  of  the  report. 
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in  the  first  annual  report.  Although  rather  extensive  literature  on  LVD 
was  known  to  exist,  it  was  recognized  at  the  outset  that  most  past  obser¬ 
vations  had  been  scattered  and  uncoordinated.  Among  the  U.S.,  European, 
and  Soviet  investigators  who  had  previously  studied  LVD  phenomena,  there 
was  little  agreement  concerning  the  validity  of  much  of  the  reported  work. 
Moreover,  there  had  been  no  very  thorough  attempts  to  resolve  the  conflicts. 
Accordingly,  the  program  began  with  an  effort  to  delineate  the  really 
important  factors  contributing  to  LVD  and  to  discard  those  concepts  and 
measurements  that  have  tended  to  confuse  the  issue.  The  objective  of 
this  initial  study  was  to  establish  the  direction  to  be  taken  in  planning 
and  executing  the  rest  of  the  investigation. 

Each  of  the  surprisingly  numerous  existing  conceptual  models  of  LVD 
initiation  and  propagation  was  first  examined  against  the  background  of 
accumulated  data.  Many  concepts  were  thus  shown  to  be  highly  improbable 
and  unworthy  of  further  study.  The  remaining  concepts,  which  seemed 
reasonably  consistent  with  experimental  information,  were  subjected  to 
fairly  simple  analyses  as  a  further  test  of  their  validity  and  relevance 
for  subsequent  research  in  this  program. 

Apart  from  quantitative  disagreements  with  available  data,  most 
models  failed  to  explain  certain  important  and  experimentally  well- 
established  physical  characteristics  of  LVD,  such  as  the  marked  influence 
of  container  geometry  and  material  (strength)  on  the  initiation  and 
stability  of  LVD  waves.  Two  conceptual  models  that  had  been  proposed 
seemed  qualitatively  consistent  with  these  behavior  characteristics  and 
therefore  deserved  careful  consideration. 

5 

Amster  et  al.  suggested  the  possibility  of  hot-spot  initiation 
induced  by  a  Mach  disk  (local  normal  shock)  occurring  at  the  center  of 
the  container  when  precursor  shocks  generated  by  container  wall  deflections 
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converge  along  the  centerline.  Although  this  model  affords  a  plausible 

explanation  of  certain  behavioral  characteristics,  a  recent  analysis  at 

Stanford  Research  Institute  (SRI)  indicates  that  the  Mach  disk  generally 

is  not  strong  enough  to  create  the  local  temperatures  needed  for  initiation. 

The  hot-spot  initiation  concept  was  retained  in  a  model  first  suggested 
6 

by  Watson  et  al.  Their  model,  which  qualitatively  fits  existing  experi¬ 
mental  observations,  envisioned  the  appearance  of  cavitation  in  the 

liquid,  followed  by  a  shock  wave  that  collapses  the  bubbles,  causing 

7 

local  hot  spots  that  act  as  initiation  centers.  Watson  had  also  taken 

photographs  showing  cavitation  ahead  of  the  reaction  zone,  and  Watson 
8 

et  al.  had  more  recently,  reported  measurements  of  elastic  precursor 

wall  waves  that  further  substantiate  the  cavitation  model  of  LVD.  Based 

4  9 

upon  the  work  done  at  SRI  under  this  program  ’  and  in  light  of  the 
properties  of  LVD  observed  to  date  by  others,  the  cavitation  model 
appeared  to  offer  the  most  plausible  mechanism  for  the  propagation  of  LVD. 

The  cavitation  model  for  LVD  assumes  that  container  interaction  with 
the  energetic  liquid  produces  cavitation  ahead  of  a  shock  in  the  liquid. 

The  liquid  shock  subsequently  compresses  the  bubbles,  causing  initiation 
and  subsequent  burning.  The  approach  taken  in  this  pi’ogram  to  develop  a 
better  theoretical  understanding  of  the  LVD  phenomenon  has  been  to  examine 
the  mechanisms  by  which  cavitation  is  produced  in  the  liquid  and  theo¬ 
retically  to  model  the  reaction  zone  in  an  attempt  to  predict  observed 
LVD  behavior.  This  theoretical  work  was  supplemented  by  experiments 
designed  to  test  the  analytical  model  critically.  In  the  following,  the 
work  on  the  cavitation  production  mechanism,  the  theoretical  modeling  of 
LVD,  and  the  experimental  results  are  summarized. 


Ill  LVD  COMPUTER  MODEL  -  WALL  INTERACTION  PRODUCING  CAVITATION 

The  validity  of  the  cavitation  model  of  LVD  requires  that  a  mechanism 

exists  for  the  production  of  cavitation  in  the  liquid  ahead  of  the  LVD 

wave.  As  a  part  of  the  current  program,  an  analysis  of  the  liquid-wall 

interaction  was  performed  to  determine  whether  significant  cavitation 

should  indeed  occur.  In  a  simplified  preliminary  analysis,  it  was  found 

that  precursor  shocks  in  the  container  walls  can  cause  a  container  wall- 

liquid  interaction  capable  of  producing  a  substantial  degree  of  cavitation 
4 

in  the  liquid.  This  analysis  was  carried  further  by  a  computer  calcu¬ 
lation  using  SRI's  two-dimensional  shock  code  to  calculate  the  stresses 
produced  in  a  liquid  contained  in  a  cylindrical  steel  pipe  or  tube  when 
shocked  at  one  end.  Large  tensile  stresses  were  produced  in  the  liquid 
across  its  entire  cross  section  by  the  wall  interaction.  The  results 
show  that  large  degrees  of  cavitation  can  be  produced  in  the  liquid. 

This  cavitation  and  the  resulting  possibility  of  hot-spot  initiations 
are  very  important,  not  only  for  LVD,  but  for  the  initiation  of  HVD  in 
liquid  (and  even  some  solid)  explosives. 

The  calculations  of  stress  wave  propagation  down  a  liquid-filled 

cylindrical  tube,  impacted  at  one  end,  were  made  using  SRI's  FIBROUS 

computer  code.  FIBROUS  is  a  two-dimensional  finite-difference  code 

10 

modeled  after  the  description  of  Wilkins,  that  can  handle  both  elastic- 
plastic  and  purely  hydrodynamic  constitutive  relations. 

To  study  further  the  cavitation  model  of  LVD,  the  stress  history  was 
computed  for  a  liquid  explosive  contained  in  a  cylindrical  steel  container 
subsequent  to  its  being  impacted  at  one  end.  The  results  of  the  compu¬ 
tations  show  that  coupling  between  the  container  wall  and  the  liquid 
produces  precursor  tension  zones  in  the  liquid  leading  to  cavitation 
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ahead  of  the  main  liquid  shock  and  that  recompression  of  the  cavitated 
zones  by  the  liquid  shock  is  sufficient  to  initiate  reaction  to  support 
the  propagation  of  the  liquid  shock.  The  simulation  for  ethyl  nitrate 
in  a  1-inch  I.D.  by  1/8-inch  wall  thickness  steel  pipe  resulted  in 
detonation  velocities  of  a^put  1.8,  2.4,  and  3.0  mm/psec,  respectively, 
for  energy  releases  of  0,  „  0.4,  and  1.0  times  the  heat  of  reaction. 

Since  experimental  LVD  propagates  at  about  2  mm/psec,  this  was  taken  as 
evidence  of  partial  energy  release  in  LVD.  The  results  of  this  work 
showed  the  details  of  the  cavitation  production  and  recompression  and 
further  substantiated  the  cavitation  model  of  LVD. 

The  computer  model  of  LVD  provides  a  means  for  studies  of  initiation 
.id  propagation  of  LVD  in  various  confinement  systems  and,  in  addition, 
has  provided  valuable  information  about  the  initiation  and  propagation 
of  LVD.  The  results  of  the  computer  modeling  of  LVD  are  fully  described 
in  Appendix  A. 
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IV  HYDROTHERMODYNAM ICS  OF  THE  CAVITATION  MODEL  OF  LVD 


Our  cavitation  model  for  LVD  is  a  reactive  shock  propagating  in  a 
cavitated  liquid.  Propagation  is  thus  a  very  complex  phenomenon  involving 
the  formation  of  a  cavitation  field  and  its  subsequent  behavior  under 
shock  compression.  The  most  important  problem  in  a  theoretical  treatment 
of  LVD  is  that  of  calculating  the  velocity  of  propagation.  Such  calcu¬ 
lations  must  take  account  of  liquid-container  wall  interaction,  and  the 
dynamics  and  thermodynamics  of  cavitation  bubbles  with  heat  transfer, 
mass  transfer,  and  chemical  reaction.  The  simplest  treatment  of  these 
processes  was  used  in  the  initial  work  to  develop  a  method  of  calculating 
LVD  velocity.  Our  approach  was  motivated  by  the  fact  that  the  classical 
C-J  hypothesis  is  satisfactory  for  calculating  the  detonation  velocity  in 
gases  even  though  a  gaseous  detonation  wave  consists  of  a  complex  system 
of  interacting  shock  waves.  In  the  simplest  treatment  of  LVD,  the  cavita¬ 
tion  field  in  virgin  liquid  is  formed  quickly  through  liquid-container 
interactions,  and  is  accounted  for  only  by  the  initial  conditions  for 
shock  propagation;  the  compressed  cavitation  field  is  assumed  to  be  in 
mechanical  equilibrium  and  to  form  a  steady-state  reaction  zone  terminating 
in  a  state  of  complete  reaction. 

It  was  demonstrated  that  the  LVD  velocity  is  uniquely  determined  by 
our  assumptions.  With  our  assumption  of  cavitation,  the  shock  is  repre¬ 
sented  by  a  Rayleigh  lino  passinr  through  the  initial  condition  of  the 
cavitated  liquid  and  not  the  initial  condition  of  the  liquid.  With  the 
assumption  of  mechanical  equilibrium  in  the  compressed  mixture  of  bubbles 
aid  liquid,  there  is  a  well-defined  pressure  at  each  point  of  the  reaction 
zone,  and  the  initial  state  for  the  burning  is  the  initial  shocked  con¬ 
dition  in  the  cavitated  liquid.  With  the  assumption  of  complete  reaction, 
the  end  of  the  reaction  zone  lies  on  the  equilibrium  products  Hugoniot 
curve  centered  on  the  initial  shocked  condition  of  the  cavitation  liquid. 
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For  the  initial,  simplest  treatment  of  the  cavitation  model,  LVD 
propagation  velocity  is  determined  by  the  C-J  condition  on  the  equilibrium 
products  Hugoniot  curve  centered  on  the  cavitated  liquid. 

The  TIGER  code  was  used  to  calculate  detonation  parameters  in  ethyl 
nitrate  for  different  degrees  of  cavitation,  and  the  results  of  these 
calculations  are  shown  in  Table  1.  The  internal  energy  of  the  cavitated 
liquid  ahead  of  the  shock  was  assumed  to  be  the  same  as  that  of  the 
uncavitated  liquid.  Calculations  with  the  BKW  and  virial  equations  of 
state  do  not  agree  with  experiment;  propagation  velocities  calculated 
with  the  ideal  gas  equation  of  state  agree  with  experiment,  but  detona¬ 
tion  pressures  lie  in  the  kilobar  region  where  this  equation  of  state 
is  inapplicable. 

The  results  of  the  calculations  led  to  the  conclusion  that  the  initial 
assumptions  were  too  restrictive.  A  more  sophisticated  treatment  of  the 
compressed  cavitation  field  was  required  to  model  LVD.  Thus,  a  more 
detailed  study  of  bubble  collapse,  heat  and  mass  transfer,  and  chemical 
reaction  was  undertaken  to  achieve  this  objective.  In  addition,  a 
modification  of  the  TIGER  code  to  allow  partial  energy  release  was  used 
to  examine  a  less  restrictive  model  and  to  determine  the  amount  of  reac¬ 
tion  required  to  support  LVD.  The  hydrothermodynamic  treatment  of  the 
cavitation  model  of  LVD  is  discussed  in  greater  detail  in  the  manuscript 
of  a  technical  paper  attached  as  Appendix  D. 


10 


Detonation  Parameters  in  Ethyl  Nitrate  Calculated  with  the  TIGER  Code 
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V  BUBBLE  DYNAMICS 

In  the  initial  study  of  bubble  dynamics,  the  characteristic  times 
for  bubble  collapse,  heat  transfer  within  a  bubble,  and  heat  transfer 
to  the  liquid  were  compared  as  a  means  of  assessing  the  relative  rates  of 
relaxation  processes  occurring  in  the  compressed  cavitation  field  behind 
the  shock.  The  collapse  of  a  bubble  was  then  treated  in  more  detail  in 
an  attempt  to  determine  conditions  for  the  onset  of  reaction  and  the 
initiation  of  LVD. 

The  dynamics  of  the  collapse  of  a  vapor  bubble  were  analyzed  to 
study  the  hot-spot  initiation  mechanism  essential  to  the  cavitation  model 
of  LVD.  The  temperature  distributions  in  a  collapsing  spherical  vapor 
bubble  and  in  the  surrounding  liquid  were  computed.  Gas-phase  reaction, 
heat  and  mass  transfer  at  the  bubble  wall,  and  motion  of  the  bubble  wall 
were  included  in  the  computations  in  order  to  simulate  closely  the  thermal 
history  of  a  collapsing  vapor  bubble  and  the  criteria  for  hot-spot 
initiation.  A  vapor  bubble  was  studied  since  during  the  very  rapid 
cavitation  that  occurs  in  liquids  undergoing  LVD,  even  if  dissolved 
gases  such  as  air  are  present,  there  is  not  sufficient  time  for  diffusion 
into  the  resulting  cavity.  Thus,  cavitated  zones  under  such  conditions 
would  consist  of  vapor  bubbles  jontaining  negligible  amounts  of  gas. 

This  work  on  bubble  collapse  loading  to  hot-spot  initiation  is 
discussed  in  detail  in  Appendix  D. 
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VI  EXPERIMENTAL  RESULTS 


As  a  part  of  the  studies  of  LVD  undertaken  in  this  program,  three 
separate  sets  of  experiments  were  carried  out  to  substantiate  the  theo¬ 
retical  model  and  to  further  delineate  the  important  mechanisms  that 
occur  in  the  initiation  and  propagation  of  LVD.  The  objective  of  the 
first  experiment  was  to  study  the  role  of  confinement,  that  of  the  second 
set  of  experiments  was  to  measure  the  wall  stress  Just  ahead  of  the 
detonation  wave,  and  in  the  third  set  of  experiments  detailed  measure¬ 
ments  were  made  of  the  stress  in  the  liquid  ahead  of  and  during  the 
passage  of  the  leading  LVD-associated  shock  wave.  The  detailed  descrip¬ 
tion  of  these  experiments  and  the  results  obtained  are  described  in  the 
first  and  second  annual  reports  and  Appendix  E. 

In  summary,  the  first  set  of  experiments  wore  run  using  metriol 
trinitrate  (TMETN)  in  steel  tubes  of  various  diameter  and  wall  thicknesses. 
For  each  configuration  LVD  velocity  was  measured  down  the  tube  for  various 
input  shock  strengths.  Once  the  proper  initiating  shock  was  determined 
for  a  nominal  1/8-inch  wall  tube,  the  wall  thickness  was  varied  keeping 
all  other  conditions  constant.  Weaker  confinement  generally  resulted  in 
decay  of  the  LVD  wave  as  it  propagated  down  the  tube,  and  stronger  con¬ 
finement  resulted  in  transition  to  UVD.  The  important  role  of  circum¬ 
ferential  or  hoop  stress  in  the  cavitation  mechanism,  as 
predicted  by  theory,  was  confirmed  by  experiments  using  a  tube  with  an 
axial  slit  down  one  side. 

In  the  second  set  of  experiments,  a  limited  number  of  shots  wore 
made,  again  using  TMETN,  to  measure  wall  stress  ahead  of  the  main  LVD 
shock.  Deformation  and  ringing  of  the  wall  were  observed,  also  in 
agreement  with  theory. 
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The  most  recent  experiments,  designed  to  measure  the  pressure  in 
the  liquid  ahead  of  and  during  passage  of  the  LVD  wave,  were  run  using 
several  liquids  and  with  1-inch  and  3/4-inch  I.D.  steel  tubes.  The 
measured  pressures  were  in  close  agreement  with  the  pressures  calculated 
by  the  computer  model.  Appendix  E  provides  a  detailed  discussion  of  this 
work. 

In  general,  the  close  agreement  between  experimental  results  and 
theoretical  predictions  based  upon  the  cavitation  model  of  LVD  further 
supports  the  cavitation  model  and  provides  a  better  understanding  of  the 
LVD  phenomenon.  In  addition,  this  capability  gives  a  good  basis  for 
further  research. 
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ABSTRACT 


To  study  the  cavitation  nodel  of  low-velocity  detonation  (LVD) ,  the 
stress  history  is  computed  for  a  liquid  explosive  contained  in  a 
cylindrical  steel  container,  subsequent  to  its  being  impacted  at  one 
end.  It  is  shown  that  coupling  between  the  container  wall  and  the 
liquid  produces  precursor  tension  zones  in  the  liquid  leading  to  cavi¬ 
tation  ahead  of  the  main  liquid  shock.  Recompression  of  the  cavitated 
zones  by  the  liquid  shock,  in  accordance  with  the  cavitation  LVD  model, 
is  assumed  to  cause  hot-spot  initiation  of  the  liquid,  and  the  sub¬ 
sequent  energy  release  in  turn  supports  the  liquid  shock.  The  simula¬ 
tion  for  ethyl  nitrate  in  a  1-inch  I.D.  by  1/8-inch  wall  thickness  steel 
pipe  resulted  in  detonation  velocities  of  about  1.8,  2.4,  and  3.0  mm/jxsec, 
respectively,  for  energy  releases  of  0.1,  0.4,  and  1.0  times  the  heat  of 
reaction.  Since  experimental  LVD  propagates  at  about  2  mm/p.sec,  this  is 
taken  as  evidence  of  partial  energy  release  in  LVD.  The  results  show 
the  details  of  cavitation  and  recompression  and  further  substantiate 


the  cavitation  model  of  LVD. 


Introduction  and  Background 


Low-velocity  detonation  (LVD)  is  a  low-order  detonation  phenomenon 
that  is  observed  in  liquid  explosives  under  particular  conditions  of 
confinement.  In  1919  LVD  was  first  observed  in  nitroglycerine  by 
Stettbacker. 1  Until  the  past  decade,  however,  the  phenomenon  remained 
relatively  obscure  and  there  was  little  agreement  among  investigators 
concerning  the  observed  properties  of  LVD.  More  recently,  because  of 
the  increasing  uses  of  highly  energetic  liquids  as  propellants,  fuels, 
and  explosives— many  of  which  are  known  to  be  susceptible  to  LVD— a 
number  of  investigations  of  LV  D  have  been  undertaken.  As  a  result  of 
these  studies,  the  important  characteristics  of  LVD  are  becoming  well 
defined,  and  plausible  explanations  of  the  observed  phenomena  have 
emerged. 

Experiments  indicate  that  LVD  can  occur  in  almost  any  highly 
energetic  monopropellant.  Typically,  LVD  propagates  at  speeds  only 
slightly  greater  than  the  liquid  sound  speed  or  in  the  range  of  1.5  to 
2.5  mm/|isec,  Stablj  LV3  is  in  general  observed  only  in  strongly  confined 
explosive  systems  and  only  i.'  the  sound  speed  of  the  container  wall 
material  exceeds  that  at  the  undisturbed  liquid.  On  the  other  hand, 
in  the  very  weak  containers  and  in  containers  with  wall  sound  speeds 
less  than  that  of  <he  liquid,  an  unstable  LVD  is  sometimes  observed 
which  propagates  in  a  pulsating  manner  with  an  average  speed  that  is 
subsonic  relative  t >  the  liquid.  Unstable  LVD  typically  propagates 
only  a  short  distance  before  dying  out, 

Though  both  tyt>es  of  detonation  are  observed  in  some  liquids,  the 
observed  properties  of  LVD  are  distinctly  different  from  those  of  HVD. 
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In  addition,  energetic  liquids  are  .generally  much  mere  sensitive  to  LVD 

initiation;  LVD  may  be  initiated  by  shock  overpressures  on  the  order  of 

1  to  5  kbar  while  shocks  of  50  to  100  kbar  are  usually  required  to 

initiate  HVD.  Moreover,  above  a  certain  minimum  charge  diameter,  HVD 

propagation  characteristics  are  independent  of  confinement  whereas  LVD 

is  strongly  dependent  upon  both  confinement  geometry  and  tin.  properties 

of  the  wall  material.  Thus,  LVD  represents  a  stable,  low-order  detonation 

wave,  the  HVD  wave.  Over  the  years,  a  number  of  models  of  the  LVD 

phenomenon  have  been  proposed  but,  until  recently,  attempts  to  explain 

1 

LVD  theoretically  have  met  with  little  success.  Stettbacker,  in  an 

attempt  to  explain  his  early  observations,  and  later  Deserhkovich  and 
2 

Andreev  suggested  tnat  the  two  detonation  velocities  in  nigroglycerine 

were  related  to  two  known  crystalline  forms  of  the  solid.  More  recently, 

3 

in  a  theoretical  paper,  Bolkhovitinov  attempted  to  explain  LVD  on  the 

basis  that  a  phase  transition  from  liquid  to  solid  takes  place  behind 

'he  initiating  shock.  In  addition,  as  a  part  of  their  study  of  hot-spot 

4 

initiation  in  explosives,  Bowden  and  Gurton  suggested  that  energy  losses 
due  to  lateral  expansion  in  the  reaction  zone  give  rise  to  LVD  behavior. 
These  theories,  although  predicting  a  low-order  detonation,  do  not  account 
for  the  experimentally  observed  characteristics  of  LVD,  particularly 
that  the  sound  speed  in  the  container  material  must  be  greater  than 
that  in  the  liquid. 

In  contrast  to  the  above  phenomenological  theories,  a  model  of  LVD 

5 

based  on  purely  thermodynamic  theory  by  Eyring  et  al.  and  a  similar 

6 

approach  taken  later  b^  Evans  utilized  a  var.lable-roaction-zone  length 
which  depended  upon  detonation  velocity.  These  models  predict  the 
existence  of  two  detonation  velocities  but  the  resulting  LVD  velocity 
varies  inversely  with  charge  diameter,  a  fact  which  is  not  supported  by 
experimental  evidence.  An  attempt  to  explain  the  LVD  phenomenon  with  a 
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semiempirical  model  based  on  varying  energy  release  rates  behind  the 

7 

initiating  shock  was  given  by  Schall  who  suggested  that  disturbed  gas 
bubbles  distributed  in  the  liquid  might  account  for  the  energy  release 
rates  responsible  for  LVT.  While  each  thermodynamic  model  of  LVD  predicts 
a  low-order  detonation  velocity,  they  all  fail  to  explain  the  experi¬ 
mentally  observed  characteristics  of  LVD,  especially  that  LVD  is  observed 
in  many  liquid  compounds  with  different  physical  properties,  propagates 
with  approximately  the  same  detonation  velocity  independently  of  the 
compound,  and  is  apparently  strongly  coupled  to  the  confinement  system. 

Recognition  of  the  close  interrelationship  between  LVD  properties, 

such  as  the  ease  of  initiation,  stability,  and  detonation  speed,  and  the 

charge  confinement  characteristics,  such  as  container  geometry,  strength, 

anil  elastic  properties  of  the  wall  material,  has  led  to  models  of  LVD 

which  include  detonation  wave-wall  coupling.  One  model,  for  example, 

that  accounted  for  the  empirical  interrelationship  between  LVD  incidence 
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and  container  wall  characteristics  was  given  by  Woolfolk  and  Amster  who 
suggested  that  the  container  wall  shock  causes  pressure  waves  in  the 
liquid  which  converge  at  the  center  to  produce  a  Mach  disk  capable  of 
initiating  detonation.  However,  in  this  model,  the  wall  shock  and 
therefore  the  Mach  disk  would  continuously  outdistance  the  LVD  wave 
unless  the  LVD  were  moving  at  a  speed  higher  than  that  observed. 

Another  empirical  model  which  also  utilizes  coupling  between  the 

wall  and  the  detonation  and  which  scorns  tc  fit  the  experimentally  observed 
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facts  was  recently  suggested  by  Watson  et  al.  and  independently  by 

10 

Voskoboinikov  et.  al  .  In  this  model,  precursor  waves  in  the  container 
wall  cavitate  the  liquid  ahead  of  a  liquid  shock.  The  cavitation  bubbles 
are  subsequently  compressed  by  the  liquid  shock  causing  a  hot-spot- 
initiated  reaction,  tnd  the  resulting  energy  release  in  turn  drives  the 
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shock.  This  model  has  been  further  verified  by  Watson  who  has  taken 

photographs  of  cavitation  ahead  of  the  reaction  zone  and  by  Watson 
12 

et  al.  who  more  recently  have  reported  tentative  measurements  of 
precursor  wall  waves  capable  of  producing  cavitation  ahead  of  the  reaction 
zone.  In  addition,  the  existence  of  individual  burning  sites  rather 
than  a  continuous  l’eaction  front  has  been  observed  by  Watson*1  in  X-ray 
photographs  of  the  LVD  inaction  wave,  and  Gibson  et  al.1:*'*5  have 
experimentally  verified  the  dependence  of  bubble  ignition  on  the  strength 
of  the  initiating  shock.  The  cavitation  model  also  explains  the  rela¬ 
tively  long  dark  region  observed  between  the  shock  front  and  the  lumi- 
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nescent  reaction  zone.  Thus,  the  cavitation  model,  though  it  involves 
an  incompletely  understood  interaction  between  the  container  and  the 
liquid  explosive  to  produce  cavitation  and  subsequent  hot-spot  reaction 
and  as  yet  unknown  criteria  for  stability,  appears  to  offer  the  most 
plausible  mechanism  for  the  propagation  of  LVI). 

Our  purpose  here  is,  on  the  basis  of  the  cavitation  model,  to 
examine  cavitation  and  recompression  mechanisms  in  both  the  initiation 
an,,  propagation  of  LVD.  A  computer  simulation  of  the  cavitation  model 
of  LVD  is  utilized  to  study  the  stress  history  in  a  system  consisting 
of  a  liquid  explosive  in  a  cylindrical  container  subsequent  to  being 
shocked  at  one  end.  The  hot-spot-initiated  reaction  is  simulated  by 
energy  release  when  a  prccavitatcd  zone  in  the  liquid  is  recompressed 
above  a  threshold  level  by  the  main  liquid  shock.  The  results  help 
clarify  the  liquid-wall  interaction  mechanism  that  leads  to  precursor 
cavitation  and  the  subsequent  collapse  of  the  cavities  by  the  liquid 
shock,  and  further  substantiates  the  cavitation  model  of  LVD. 


Cavitation  Model  of  LVD 


Experimental  sutdies  of  LVD  are  usually  carried  out  in  an  arrange- 
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mcnt  similar  to  that  used  to  study  the  sensitivity  of  liquid  explosives. 
The  liquid  to  he  studied  's  placed  in  a  container,  most  often  a  cylindri¬ 
cal  tube,  which  is  shocked  at  one  end  by  an  explosive  donor.  The  strength 
of  the  initiating  shock  can  be  varied  by  using  attenuators  (usually 
Lucite)  of  different  thicknesses.  For  a  number  of  liquids,  initiating 
shocks  on  the  order  of  1  to  10  kbar  typically  produce  LVD  whereas  much 
stronger  shocks,  on  the  order  of  100  kbar,  directly  initiate  UVD. 
Experimenters  in  the  past  have  studied  LVD  phenomena  with  high-speed 

cameras,  velocity-measuring  probes,  and  strain  gages  to  measure 
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circumferential  wall  strain. 

When  the  initiating  donor  shock  wave  impacts  the  explosive  system, 
the  shocks  produced  propagate  down  the  cylinder  from  the  donor  end.  If 
the  sound  speed  in  the  container  wall  material  is  greater  than  that  in 
the  liquid,  then  the  wall  shock  outdistances  the  liquid  shock  and  causes 
a  precursor  wave  system  ahead  of  the  liquid  shock.  Calculations  discussed 
in  the  following  sections  show  that  those  precursor  waves  cause  zones  of 
tension  capable  of  cavitating  the  liquid.  The  cavitated  liquid  is  sub¬ 
sequently  recompressed  by  the  liquid  shock  leading  to  the  possibility  of 
hot-spot  initiation.  This  mechanism  is  probably  the  manner  in  which  LVD 
initiation  takes  place  and  under  some  conditions  as,  for  example,  for  a 
stronger  donor  shock,  could  lead  to  UVD  initiation.  The  calculations 
demonstrate  that  relatively  weak  donor  shocks,  on  the  order  of  1  kbar, 
are  sufficient  to  produce  significant  cavitation  in  the  liquid,  thus 
accounting  for  the  high  shock  sensitivity  of  liquid  explosives  to  hot¬ 
spot  initiation  and  to  LVD. 


The  cavitation  model  of  LVD  requires  that  precursor  wall  waves 
produce  tension  zones  leading  to  cavitation  ahead  of  the  liquid  shock. 

If  the  wall  shock  caused  by  the  donor  were  the  only  mechanism  for  cavi- 
tating  the  liquid,  then  all  LVD  waves  would  be  inherently  unsteady,  since 
in  the  experimental  conditions  where  stable  LVD  is  observed,  the  wall 
shock  continually  outdistances  the  LVD  waves.  However,  if  the  wall 
deflection  due  to  the  LVD  wave  itself  is  capable  of  producing  precursor 
cavitation,  then  a  self-sustaining,  steady  LVD  is  possible.  Though 
Watson  et  al,^*  report  some  preliminary  experimental  evidence  that  there 
is  a  precursor  pressure  wave  attached  to  the  reaction  zone,  it  has  not 
been  established  that  such  a  mechanism  leads  to  precursor  cavitation. 

In  either  case,  it  is  clear  that  when  the  sound  speed  in  the  container 
wall  material  is  less  than  that  in  the  liquid,  precursor  wall  waves 
caused  by  either  the  donor  shock  or  the  LVD  wave  are  not  possible. 

Pulsating  LVD  is  sometimes  observed  when  the  wall  sound  speed  is 
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only  slightly  greater  than  that  of  the  liquid;  it  may  occur  when  the 
complex  interaction  between  the  donor  and  the  receptor  causes  hot-spot 
initiation  which  starts  to  accelerate  toward  a  full  LVD,  Subsequently, 
as  the  resulting  reaction  wave  approaches  the  .all  sound  speed  there  is 
no  mechanism  for  the  generation  of  precursor  cavitation  and  the  reaction 
dies  out  from  the  lack  of  initiation  centers.  As  the  reaction  zone  slows 
down  due  to  decreased  energy  release  rate,  the  wall  wave  may  again  out¬ 
distance  it,  once  again  causing  cavitation  ahead  of  the  reaction  zone. 

The  whole  process  may  then  repeat  itself  in  an  oscillatory  fashion 
resulting  in  pulsatory  LVD.  Clearly,  whether  such  an  unstable  LVD  will 
continue  to  propagate  must  depend  critically  upon  the  experimental  con¬ 
ditions.  In  fact,  unsteady  LVD  is  frequently  observed  to  die  out 
altogether. 
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Computer  Simulation 


A  computer  study  of  cavitation  and  recompression  in  an  energetic 
liquid  contained  in  a  cylindrically  symmetric  tube  impacted  at  one  end 
was  undertaken  in  an  attempt  to  simulate  the  initiation  and  propagation 
of  LVD.  Ethyl  nitrate,  contained  in  a  1-inch  I.D.  by  1/8-inch  wall 
cylindrical  steel  tube,  was  chosen  for  the  calculations  as  being  rep¬ 
resentative  of  a  typical  experimental  setup  used  to  study  the  shock 
initiation  of  LVD  in  the  practical  part  of  the  research  program.  Of 
particular  interest  in  these  calculations  are  the  following  factors: 

(1)  Does  the  stress  wave  traveling  down  the  pipe  walls  create  a 
region  of  tension  in  the  liquid  sufficient  to  produce 
cavitation  ahead  of  the  main  compression  wave  in  the  liquid? 

(2)  What  is  the  propagation  velocity  of  LVD  and  how  does  it 
depend  upon  the  amount  of  energy  released  in  the  reaction 
zone? 

(3)  Is  a  stable  LVD  wave  (one  that  propagates  at  a  constant 
velocity)  established  at  some  distance  (equal  to  several 
pipe  diameters)  down  the  pipe? 

Computer  simulations  wore  performed  using  the  SRI  FIBROUS  code. 

This  code  is  a  two-dimensional  finite  difference  stress  wave  propaga- 
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tion  code  modeled  after  the  description  of  Wilkins,  that  includes 
both  elastic-plastic  and  hydrodynamic  constitutive  relations.  When  the 
FIBROUS  code  is  used  in  the  axially  symmetrical  geometry,  a  rectangular 
grid  parallel  to  the  axis  of  the  cylinder  extending  from  the  axis  to 
the  edge  of  the  cylinder  is  specified,  as  shown  in  Figure  1.  The  code 
treats  each  cell- of  this  plane  grid  as  if  it  contained  the  volume  swept 
out  by  the  360-degree  rotation  of  the  grid  about  the  cylindrical  axis. 
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A  similar  projectile  cylinder  impacts  the  target  from  the  right  to  begin 
the  computation.  (The  shock  produced  in  the  explosive  system  by  the 
projectile  cylinder  simulates  the  shock  delivered  in  the  experimental 
setup  by  the  donor-attenuator.)  At  successive  time  increments  after 
impact,  all  of  the  relevant  parameters  needed  to  characterize  the  stress 
wave  propagation  are  calculated  for  each  grid  cell.  These  include 
position,  velocity,  triaxial  stresses,  hydrodynamic  pressure,  internal 
energy,  and  specific  volume.  Certain  specified  parameters  at  designated 
cells  can  be  stored  by  the  computer  at  each  time  increment,  and  these 
can  be  tabulated  or  plotted  at  the  end  of  the  calculation  to  give,  for 
example,  stress  and  specific  volume  histories,  or  stress-volume  loading 
and  unloading  paths. 

For  each  cell  during  each  cycle  the  calculation  proceeds  basically 

as  follows:  the  stresses  in  the  adjacent  cells  during  the  previous  cycle 

are  used  to  compute  the  current  acceleration  for  the  cell  which,  along 

with  its  previous  velocity  and  position,  is  then  used  to  compute  its 

current  velocity  and  position.  From  the  latter  are  obtained  the  current 

triaxial  strains,  distortional  strain  energy,  and  the  specific  volume 

which— together  with  the  equation  of  state  or  constitutive  relation  of 

the  material  in  the  cell--are  used  to  compute  the  current  triaxial 

stresses,  the  hydrostatic  pressure,  and  the  internal  energy.  The 
»i  IV 

Mie-Gruneisen  equation  of  state  is  used  to  relate  changes  in  the  internal 
energy  to  changes  in  the  hydrostatic  pressure.  The  time  increment  between 
successive  cycles  is  chosen  so  that  a  stress  wave  propagates  less  than 
the  distance  between  two  cells  during  that  time. 

The  computational  grid  for  the  simulation  of  LVD  was  set  up  as  in 
Figure  1.*  The  target  grid  contained  d  rows  of  liquid  cells  and  one  row 
of  steel  pipe  cells  in  the  radial  direction;  each  row  contained  97  cells 

*  The  figures  are  shown  at  the  end  of  this  appendix. 


in  the  axial  direction.  However,  only  about  the  first  half  of  the  target 
in  t ho  axial  direction  yielded  useful  stress  histories  in  the  liquid 
cells,  since  reflection  of  the  stress  wave  in  the  pipe  wall  from  the 
free  surface  disrupted  the  stress  flow  in  the  latter  half  of  the  target. 
Therefore  the  usable  LVD  simulation  calculations  were  carried  out  only 
to  about  6  diameters  down  the  tube.  The  projectile,  which  contained  5 
rows  of  cells  in  the  radial  direction  with  each  row  containing  10  cells 
in  the  axial  direction,  impacted  the  target  at  a  velocity  of  0.5  mm/psec. 
The  cells  for  both  the  target  and  projectile  grids  were  3  mm  square. 
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An  elastic-plastic,  equation-of-state  formulation  for  Armco  iron 
was  used  for  the  constitutive  relations  for  the  steel  pipe  cells.  It 
included  a  nonlinear  hydrostatic  pressure-specific  volume  (P-V)  llugoniot 
curve  given  by 

P  =  1,590  p  +  5,170  p2  +  51,700  p‘J  , 

Vn  3 

where  p  =  —  -  1,  P  is  in  kbar  and  V  =  0.1274  cm/gm.  A  shear 

strength  of  7.9  kbar  and  a  shear  modulus  of  820  kbar  were  used.  The 

equation  of  state  for  water  was  used  for  the  liquid  cells,  inasmuch  as 

the  density  and  shock  properties  of  water  are  similar  to  those  of  the 

liquid  propellants  of  interest  in  LVD  studies.  The  hydrostatic  pressure 

of  the  liquid  was  defined  by 

2  3 

P  =  25. 6p  +  01. 2p  +  122. 7p  . 

It  was  found  that  when  a  purely  hydrodynamic  equation  of  state  (one  in 
which  the  shear  strength  is  zero)  is  used  for  the  liquid,  the  liquid 
colls  next  to  the  wall  become  highly  elongated  because  liquid  in  the 
center  of  the  tube  tends  to  flow  down  the  tube  whereas  the  liquid 
adjacent  to  the  wall,  is  prevented  from  sliding  along  the  inside  of  the 
wall  (no  sliding  interfaces  are  allowed  in  the  FIBROUS  code) .  Since  a 


highly  elongated  cell,  particularly  one  in  which  one  diagonal  is  much 
larger  than  the  other,  may  lead  to  calculational  inaccuracies,  the  liquid 
was  given  a  small  shear  strength  of  0.1  kbar.  It  is  not  expected  that 
this  had  any  significant  effect  upon  the  simulation.  The  equation  of  state 
of  C-7  was  used  for  the  constitutive  relations  for  the  projectile 
material.  It  will  not  be  described  here,  since  the  projectile  equation 
of  state  has  no  effect  upon  the  stress  history  for  time  intervals 
relevant  to  the  LVD  simulation. 

The  cavitation  of  the  liquid  ahead  of  the  liquid  stress  wave  and 
the  subsequent  detonation  of  the  cavitated  liquid  upon  arrival  of  the 
liquid  shock  were  simulated  in  the  following  manner:  whenever  the  ten¬ 
sile  stress  (pressure)  in  a  cell  of  the  liquid  exceeded  0.02  kbar,  this 
cell  was  tagged  and  the  tensile  stress  was  reset  to  0.02  kbar.  So  when 
further  tensile  strain  was  applied  to  this  cell,  the  cell  continued  to 
expand  at  a  constant  negative  pressure  level  of  0.02  kbar.  This  approx¬ 
imation  to  the  liquid  stress-strain  behavior  during  cavita¬ 
tion  was  handled  with  only  minor  adjustment  to  the  FIBROUS  code,  A 
more  realistic  treatment  of  cavitation  that  is  possible  but  involves 
more  substantial  changes  in  the  FIBROUS  code  will  be  discussed  below. 

All  of  these  tagged  cells  (which  had  undergone  negative  stress  in 
excess  of  0.02  kbar  and  hence  cavitation)  subsequently  experienced 
recompression  duo  to  the  arrival  of  the  main  compressional  wave  in  the 
liquid.  When  the  pressure  in  the  tagged  cells  exceeded  0,4  kbar  in 
compression,  an  amount  of  energy  equal  to  the  energy  that  is  expected 
to  be  released  by  the  hot-spot  reaction  was  added  to  the  internal  energy 
of  that  cell  to  simulate  the  energy  release  in  the  LVD  reaction  zone. 

The  energy  was  added  in  three  equal  amounts  over  three  calculational 
cycles  (the  time  between  successive  cycles  in  the  computation  was  about 


0.35  psec)  since  a  detonation  due  to  hot-spot  initiation  of  collapsing 

bubbles  is  not  expected  to  occur  instantaneously  throughout  one  region 

of  the  liquid.  Since  there  is  uncertainty  as  to  the  amount  of  energy 

actually  released  (specifically  there  is  evidence  that,  only  a  fraction 

of  the  heat  of  reaction  is  released  during  LVD)  ,  calculations  were  made 

using  three  different  amounts  of  energy  (EDET).  to  be  added  to  each  cell, 

in  order  to  assess  the  effect  of  fractional  energy  release  upon  the 

overpressure  and  propagation  speed  of  the  LVD  wave.  The  heat  of  reac- 
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tion  of  ethyl  nitrate  which  is  approximately  5  x  10  ergs/gin,  was  used 
as  the  basis  of  EDET.  A  recompression  of  cavitated  zones  to  the  value 
of  0.4  kbar  was  chosen  as  the  compression  threshold  for  the  initiation 
of  reaction  so  that  the  low-amplitude  compressional  waves  due  to  ringing 
behind  the  initial  wall  shock  would  not  predetonate  the  liquid.  This 
feature,  in  connection  with  the  stress  history  in  the  liquid,  is 
discussed  in  more  detail  below. 

Results  of  LVD  Computations 

The  important  events  that  occur  after  impact  are  shown  in  Figure  2 

g 

for  three  simulations  in  which  the  EDET  equaled  5  x  10  ergs/gm, 
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2  x  10s  ergs/gm,  and  5  x  10  ergs/gm  or,  respectively  approximately 
0.1,  0.4,  and  1.0  times  the  heat  of  reaction  of  ethyl  nitrate.  The 
initial  stress  wave  in  the  pipe  wall,  which  travels  at  approximately 
5.3  mm/psec,  quickly  outdistances  that  in  the  liquid  which  travels  at 
approximately  1.8  mm/psec .  When  the  former  sufficiently  outdistances 
the  latter,  or  at  approximately  0.7  I.D.  down  the  pipe,  the  following 
sequence  of  events  takes  place.  The  compi’essivc  triaxial  stresses  in 
the  pipe  wall  cause  a  smn.ll  compressional  wave  to  propagate  into  the 
liquid  toward  the  center  of  the  pipe  and  at  the  same  time  cause  the  pipe 
wail  to  expand  radially.  With  this  increase  in  radial  strain,  the  hoop 


stress  in  the  pipe  wall  goes  into  tension  and  a  tensile  pulse  is  sent 

into  the  liquid  shortly  behind  the  initial  compressive  pulse.  As  this 

tensile  pulse  proceeds  toward  the  center  of  the  pipe,  the  liquid  begins 

to  cavitate,  first  at  the  cell  near  the  pipe  wall  and  finally  at  the 

central  cell.  Eventually  the  main  compressive  stress  wave  in  the  liquid 

reaches  the  cavitated  zones  and  recompresses  the  liquid  to  above  0.4  kbar, 

which  triggers  the  detonation.  Since  cavitation  occurs  earlier  in  the 

liquid  cells  near  the  edge  of  the  pipe  than  in  the  cells  at  the  center, 

and  since  the  main  compressional  wave  in  the  liquid  reaches  all  of  the 

cells  at  one  axial  position  at  about  the  same  time,  detonation  begins 

earliest  in  the  cell  adjacent  to  the  wall  (at  about  0.7  I.D.  down  the 

pipe)  and  latest  at  the  center  cell  (at  about  1.2  I.D.).  This  result 

is  consistent  with  the  observations  of  LVD  in  nitroglycerine  reported 
19 

by  Sosnova  et  al.  that  the  most  intense  reaction  is  initiated  in  the 
nitroglycerine  adjacent  to  the  container  wall  and  that  the  reaction  wave 
profile  is  concave  with  respect  to  the  direction  of  propagation  in  the 
initiating  region.  It  should  be  noted,  however,  that  since  the  inward- 
moving  tensile  pulses  converge  along  the  center  line,  the  tensile  stress 
is  highest  there  and  therefore  the  liquid  in  the  center  might  be  more 
highly  cavitated.  It  is  thus  possible  that  under  some  conditions, 
perhaps  when  the  initiating  shock  is  relatively  weak,  LVD  might  occur 
primarily  along  the  central  region. 

Figures  3a  and  3b  depict  the  computed  stress  histories  at  the 
central  cell  and  at  the  coll  adjacent  to  the  pipe  wall,  respectively, 
for  various  distances  down  the  tube  in  the  case  for  which  EDET  =  5  x  108 
ergs/gm.  For  each  stress  history  there  is  initially  a  relatively  small 
compressive  pulse  followed  by  a  region  of  tension  where  the  tensile 
stress  is  about  0.02  kbar.  These  wall-induced  precursor  waves  are 
followed  eventually  by  the  main  liquid  compressive  wave  which  triggers 
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the  detonation.  Note  that  for  distances  down  the  pipe  greater  than  about 
2.5  diameters,  there  are  one  or  more  additional  cycles  of  precursor- 
induced  compression  and  tension  in  the  liquid  prior  to  the  arrival  of 
the  main  liquid  wave.  These  waves  are  due  to  the  further  ringing  of  the 
pipe  following  the  initial  compressional  wave.  If  the  compressional 
pulses  following  the  initial  cavitation,  but  prior  to  the  arrival  of 
the  main  liquid  compression  wave,  were  strong  enough  to  collapse  the 
voids  and  cause  detonation,  then  the  LVD  wave  would  be  coupled  directly 
to  the  pipe  wall  wave  and  would  travel  at  approximately  the  wave  velocity 
in  the  pipe  rather  than  at  a  velocity  nearly  that  of  the  liquid  shock. 
Since  it  is  known  experimentally  that  this  does  not  happen,  detonation 
was  not  triggered  until  the  main  liquid  shock  arrived.  A  detonation 
threshold  of  0.4  kbar  was  chosen,  which  in  the  cases  investigated  is 
higher  than  the  magnitude  of  the  pressure  in  the  liquid  caused  by 
ringing  of  the  pipe  wall. 

Figures  4a  and  4b  depict,  respectively,  the  computed  axial  and 
circumferential  (hoop)  stresses  in  the  pipe  at  two  distances  down  the 
tube  for  the  same  computation  as  above.  The  axial  stress  history  of 
Figure  4a  clearly  shows  the  initial  wall  shock  traveling  at  a  speed  of 
about  5.3  min/psec  followed  by  a  more  gradual  stress  release  and  ringing. 
The  hoop  stress  shown  in  Figure  4b  rings  back  and  forth  between  com¬ 
pression  and  tension  following  the  initial  pipe  wave  until  a  strong 
tension  is  produced  by  the  detonation  of  the  liquid  inside  the  pipe. 

Figure  5  shows  the  stress-strain  loading  and  unloading  curves  of 
the  liquid  cell  located  at  the  center  of  the  pipe  1.2  diameters  down 
the  pipe  from  the  impacted  end.  The  effects  of  the  cavitation  and 
detonation  parameters  used  in  the  computations  can  be  seen  by  comparing 
the  two  cases  in  which  the  liquid  does  and  does  not  undergo  cavitation 


and  detonation. 
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The  simulations  in  which  EDET  of  2  x  10  and  5  x  10  ergs,  gm  were 
used  resulted  in  stress  histories  similar  to  those  depicted  in  Figures 
3  and  4.  The  two  principal  differences  were  that  the  LVD  propagated  at 
a  slightly  higher  velocity— 2.36  and  2.99  mm/psec,  respectively,  than 
the  1.84  mm/psec  for  the  run  using  the  lower  EDET,  and  the  pressure 
induced  in  the  liquid  by  the  LVD  was  significantly  higher.  Figure  6 
contains  stress  histories  at  the  central  cell  located  1.2  diameters 
down  the  pipe  for  all  three  simulations  using  the  different  values  of 
EDET  as  well  as  the  stress  history  for  a  simulation  in  which  no  detonation 
was  allowed  to  take  place.  The  simulation  represented  by  these  calcula¬ 
tions  closely  approximates  LVD  in  ethyl  nitrate  contained  in  a  1-inch 
I.D.  by  1/8-inch  wall  steel  pipe.  Since  measured  LVD  velocities  for 
this  system  are  on  the  order  of  1.9  mm/psec,  the  above  results  indicate 
that  partial  energy  release  is  likely  in  LVD.  This  hypothesis  is  further 

substantiated  by  thermodynamic  calculations  of  LVD  which  also  indicate 

20 

that  partial  energy  release  is  likely. 

An  additional  computer  simulation  was  made  using  lead  as  the  pipe 
wall  material.  The  stress  wave  velocity  in  the  lead  was  only  slightly 
higher  (2,1  mm/^sec)  than  that  in  the  liquid  for  the  particular  geometry 
and  pressures  of  the  simulation.  The  result  was  that  the  wall  wave 
never  quite  ran  far  enough  ahead  of  the  liquid  wave  to  cause  any  zone 
of  tension  in  the  liquid;  therefore,  no  cavitation  occurred  to  initiate 
LVD.  Computations  using  different  lead  pipe  wall  thicknesses  and 
stronger  initiating  shocks  could  conceivably  result  in  precursor  cavi¬ 
tation  and  an  oscillating  LVD  wave.  However,  further  computations  of 
this  typo  were  beyond  the  scope  of  the  current  work. 
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Model  for  Cavitating  Liquid 

The  model  for  cavitation  used  in  these  computations,  in  which  the 
liquid  simply  expands  at  constant  pressure  when  the  tensile  pressure 
exceeds  0.02  kbar  and  then  subsequently  contracts  at  constant  pressure 
until  the  strain  corresponds  to  that  at  0.02  kbar  in  tension,  is  only  a 
qualitative  approximation  of  the  cavitation  behavior  of  a  liquid  under¬ 
going  tension.  The  more  detailed  phenomena  of  cavitation  in  a  liquid 

undergoing  dynamic  tensile  failure  in  the  kilobar  stress  region  have 

21 

been  discussed  by  Erlich  et  al.  A  computer  model  of  the  dynamic 
equation  of  state  of  a  cavitating  liquid  was  developed  which  calculated 
growth  and  collapse  of  bubbles  in  a  liquid  strained  at  high  rates.  This 
model,  which  uses  a  viscous  growth  law  and  calculates  the  void  volume 
and  resulting  stress  relaxation  for  a  given  tensile  strain,  could  be 
inserted  in  the  FIBROUS  code  to  describe  more  quantitatively  the  cavi¬ 
tation  and  subsequent  bubble  collapse  for  each  liquid  cell  in  the  LVD 
simulation.  The  only  unknown  parameter  in  such  a  calculation  would  be 
the  void  nueleation  frequency  or  the  initial  void  volume,  which  can 
either  be  determined  to  sufficient  accuracy  for  a  given  material  by 
experiments  of  the  type  described  in  Reference  21  or  by  an  approximation 
based  on  data  for  other  liquids,  When  a  liquid  is  subjected  to  tensile 
strain,  rapid  stress  release  occurs  due  to  the  nueleation  and  growth  of 
c*  ation  bubbles.  Thus,  rather  than  using  the  strain  at  constant 
negative  stress  to  simulate  the  cavitation  behavior  as  was  done  here,  a 
more  realistic  model  would  include  dynamic  bubble  growth  and  the 
resulting  stress  release.  In  addition,  once  a  zone  of  the  liquid  has 
been  cavitatod,  a  modified  equation  of  state  to  include  the  dynamic 
bubble  motion  could  be  used.  Computations  using  those  improvements 
would  give  a  more  realistic  picture  of  the  cavitation  and  void  collapse 
phenomena  and  could  thus  greatly  improve  the  cavitation  model  of  LVD. 
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Discussion 


The  axisymmetric  computer  calculations  to  simulate  LVD  in  a 

cylindrical  geometry  show  clearly  the  manner  in  which  the  container 

wall-liquid  interaction  produces  zones  of  tension  leading  to  cavitation 

of  the  liquid.  Tension  is  produced  in  the  liquid  by  waves  propagating 

from  the  container  walls  into  the  liquid.  The  waves  propagate  into  the 

liquid  -from  the  container  walls  much  like  Mach  waves  and  converge  at  the 

tube  centerline.  If  one  computes  the  Mach  angle,  approximately  equal  to 

sin  1  (V^/VO  ,  where  a»  1.8  mm/psec  is  the  propagation  speed  of  the 

wall  waves  and  2:  5.3  mm/psec  is  the  wave  propagation  speed  in  the 

o 

liquid,  a  value  of  approximately  20  is  obtained  for  the  present  system. 

This  indicates  that  the  time  delay  for  precursor  waves  to  propagate  from 
the  tube  wall  to  the  center  cell  should  be  about  '3.3  psec,  a  result  that 
is  consistent  with  that  obtained  by  comparison  of  the  precursor  waves 
in  Figures  3a  and  3b.  Thus,  the  cavitation  produced  by  the  precursor 
wall  waves  is  a  complex  phenomenon  that  clearly  cannot  be  accounted  for 
by  a  one-dimensional  analysis. 

The  recompression  of  2  cavitated  zone  by  the  primary  liquid  shock 
is  assumed  to  initiate  hot-spot  reactions,  and  the  subsequent  energy 
release — simulated  by  dumping  an  amount  of  energy  EDET  into  the  recompressed 
zone--drives  the  primary  liquid  shock.  As  can  be  seen  by  a  comparison  of 
Figures  3a  and  3b,  once  the  detonation  is  fully  developed  the  liquid 
shock  is  approximately  planar.  Thus,  even  though  the  liquid  next  to 
the  wall  is  cavitated  prior  to  that  in  the  center  of  the  liquid,  rocom- 
pression  occurs  almost  simultaneously  over  the  tube  cross  section  once 
the  liquid  shock  is  more  than  1  or  2  diameters  from  the  initiation  end. 
Without  detailed  calculation  of  the  bubble  dynamics  and  appropriate 
mo<  cation  of  the  liquid  constitutive  equations,  the  detailed  effects 
of  ;he  time  history  upon  the  cavitation  field  cannot  be  fully  explained. 
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During  the  initiation  phase,  which  occurs  approximately  within  the 
first  two  diameters  from  the  donor  end  of  the  charge,  the  liquid  shock 
is  close  enough  behind  the  wall  shock  to  be  within  the  Mach  wave  struc¬ 
ture  causing  cavitation.  Thus,  reaction  takes  place  only  in  the  vicinity 

19 

of  the  wall,  the  region  where  initiation  was  observed  by  Sosnova  et  al. 
Watson  et  al.**  observed  a  reaction  zone  in  the  region  downstream  of  the 
initiation  region  where  cavitation  is  generated  across  the  entire  tube 
cross  section  before  arrival  of  the  liquid  shock. 

The  ringing  of  the  tube  wall  behind  the  initial  wall  shock  causes 

alternating  zones  of  tension  and  compression  in  the  liquid  ahead  of  the 

main  liquid  shock  (see  Figure  3).  The  period  of  the  oscillations  is 

about  16  psec  for  the  present  configuration.  Studies  of  liquid  failure, 
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or  cavitation  at  high  stress  rates,  show  that  significant  bubble  growth 
and  the  accompanying  stress  release  can  occur  in  less  than  1  psec  for 
stresses  on  the  order  of  those  observed  here.  It  is  therefore  likely 
that  the  ringing  stresses  produce  significant  bubble  motion  perhaps 
capable  of  collapsing  the  cavitation  produced  by  the  initial  tension 
wave.  Since  experimental  observations  show  that  the  reaction  zone  is 
not  directly  coupled  to  and  does  not  move  with  the  initial  wall  shock, 
as  would  occur  if  the  ringing  compressions  caused  bubble  collapse  and 
hot-spot  initiation,  the  artifice  was  adopted  here  that  a  compression 
greater  than  0,4  kbar  was  required  to  cause  initiation.  This  effectively 
permitted  hot-spot  initiation  at  a  precavitated  zone  only  upon  arrival 
of  the  main  liquid  shock.  If,  as  seems  likely,  ringing  of  the  liquid 
in  the  tube  ahead  of  the  main  liquid  shock  is  strong  enough  to  collapse 
the  bubbles  but  not  strong  enough  to  cause  hot-spot  initiation,  then  the 
main  liquid  shock  will  travel  into  alternating  regions  of  cavitated  and 
uncontacted  liquid.  This  is  a  very  possible  explanation  for  the  pulsating 
type  of  LVD,  In  which  the  detonation  travels  a  short  distance,  this  dies 
out,  then  begins  again  farther  down  the  tube,  and  so  forth. 
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That  the  recompression  due  to  wall  ringing  does  not  cause  initia¬ 
tion  may  be  due  to  several  reasons.  First,  the  threshold  stress  for 
hot-spot  initiation  by  a  given  bubble  .field  is  not  known.  It  is  there¬ 
fore  possible  that  the  compression  due  to  ringing  is  not  strong  enough 
to  cause  initiation.  On  the  other  hand,  the  equation  of  state  of  the 
cavitated  liquid  is  only  approximated  in  the  present  calculations  and  a 
more  detailed  model  of  the  cavitation  zone  might  provide  a  much  different 
stress  in  the  liquid  due  to  the  wall  ringing.  For  example,  the  recom¬ 
pression  by  the  final  ringing  wave  following  the  first  tension  zone  would 
be  propagating  into  a  dynamically  cavitating  liquid.  Dispersive  effects 
on  the  wave  propagation  caused  by  the  bubble  field  could  significantly 
alter  the  stress  pattern  from  that  computed  here.  The  inclusion  of  a 
more  detailed  equation  of  state  for  the  cavitated  liquid  which  accounts 
for  the  dynamic  bubble  motions  would  represent  a  considerable  extension 
of  the  current  calculations  and  is  currently  being  considered  as  a 
subject  for  further  study. 

As  noted  previously,  if  the  LVD  is  to  be  truly  a  steady-state 
phenomenon,  the  wall  wave  caused  by  the  detonation  itself  must  be 
capable  of  causing  precursor  cavitation  independently  of  the  wall  wave 
caused  by  the  initial  shock.  (A  quasi-steady  state  LVD ,  however,  could 
occur  for  a  considerable  distance  down  a  pipe  if  the  cavitation  bubbles 
remained  in  the  liquid  long  after  the  ringing  due  to  the  initial  wall 
shock  had  disappeared.)  In  the  present  calculations  there  is  no  indi¬ 
cation  that  this  is  the  case.  One  calculation  was  extended  to  an  axial 
distance  of  Z/D  =  15  and,  oven  in  the  most  energetic  case  (EDET  = 

5  x  10  ergs/gm)  in  which  the  shock  pressure  was  in  excess  of  10  kbar, 
detonation-wavo-coupled  precursor  wall  waves  capable  of  causing  cavi¬ 
tation  directly  ahead  of  the  liquid  shock  wore  not  observed.  Thus, 
while  not  conclusive,  the  work  here  indicates  that  cavitation  is  caused 
only  by  the  precursor  wall  wave  generated  by  the  initiating  shock. 
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Conclusions 


The  computations  described  above  support  the  cavitation  model  for 
the  initiation  and  propagation  of  LVD  and  provide  a  first  step  towards 
a  detailed  quantitative  description  of  these  processes.  The  results 
support  the  hypothesis  that  onl v  partial  energy  release  occurs  in  LVD. 
There  is  no  indication  that  the  wall  wave  caused  by  the  detonation  itself 
produces  zones  of  tension  capable  of  producing  cavitation.  However, 
further  conclusions  regarding  these  results  must  await  additional 
calculations  for  different  geometries,  different  materials,  and  more 
accurate  modeling  of  the  equation  of  state  of  the  cavitated  liquid. 
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FIGURE  A-1  RECTANGULAR  GRID  REPRESENTATION 
OF  LIQUID-FILLED  PIPE  USED  IN 
FIBROUS  COMPUTER  CODE 
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FIGURE  A-2  WAVE  DIAGRAM  SHOWING  RELEVANT  EVENTS 
IN  LVD  COMPUTER  SIMULATION 
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FIGURE  A-4  AXIAL  (a)  and  CIRCUMFERENTIAL 
IN  THE  PIPE  WALL  AT  SEVERAL 
SIMULATION  IN  WHICH  EDET  =  5 


CIRCUMFERENTIAL  STRESS 
HISTORIES  IN  THE  PIPE  WALL 
AT  SEVERAL  LOCATIONS  FOR 
THE  SIMULATION  IN  WHICH 
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FIGURE  A-5  STRESS-STRAIN  LOADING  AND  UNLOADING 
CURVE  FOR  THE  LIQUID  IN  THE  LVD 
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FIGURE  A-6  PRESSURE  HISTORIES  OF  THE  LIQUID  CELL 

IN  THE  CENTER  OF  THE  PIPE  AT  A  DISTANCE 
OF  1.2  INSIDE  DIAMETERS  FROM  THE  IMPACTED 
END  FOR  DIFFERENT  VALUES  OF  EDET  (EDET 
FOR  ETHYL  NITRATE  s  5  X  109  ergs/gm) 


A-27 


Appendix  B 

HYDROTHERMODYNAMICS  OF 

THE  CAVITATION  MODEL  FOR  LOW  VELOCITY  DETONATION 


i 


HYDROTHERMODYNAMICS  OF 

THE  CAVITATION  MODEL  FOR  LOW  VELOCITY  DETONATION* 

By 

M.  Cowperthwaite 
Physical  Sciences  Division 
Stanford  Research  Institute 
Menlo  Park,  California 


$ 

This  work  was  supported  by  the  U.S.  Air  Force  Office  of  Scientific 
Research  under  Contract  No.  F44620-69-C-0079. 


B-iii 


Introduction 


Shock  wave  studies  during  the  last  decade,  extending  earlier 

7 

work  discussed  in  the  book  by  Taylor,  have  established  that  many  liquid 
propellants  and  explosives  exhibit  two  modes  of  self-sustaining  detona¬ 
tion.  These  modes  have  distinct  velocities  of  propagation  and  have  been 
characterized  accordingly  as  high-velocity  detonation  (HVD)  and  low- 
velocity  detonation  (LVD) .  The  high-velocity  mode  propagates  at  approxi¬ 
mately  6  mm/psec  with  a  pressure  of  about  100  kilobars;  the  low-velocity 
mode  propagates  with  a  slightly  higher  velocity  than  ambient  soundspeed, 
at  approximately  2  mm/psec  with  a  pressure  of  about  10  kilobars.  The 

hydrodynamic  theory  of  detonation,  together  with  the  Chapman-Jouguet 
8  9 

(C-J)  ’  hypothesis,  provides  a  satisfactory  model  for  calculating  HVD 
velocities  when  the  equation  of  state  of  the  detonation  products  is 
known,  but  not  for  calculating  LVD  velocities  at  the  present  time. 

The  Cavitation  Model  for  LVD 

The  present  work  is  concerned  with  the  problem  of  formulating  a 

model  for  calculating  the  propagation  velocity  of  L\D.  It  is  based  on 

the  physical  model  proposed  by  workers  at  the  Bureau  of  Mines  after  an 

extensive  experimental  investigation  of  LVD.  This  model  was  discussed 
5 

by  Watson  at  the  Twelfth  International  Symposium  on  Combustion.  LVD 
is  modeled  as  a  reactive  shock  propagating  in  a  cavitated  liquid.  The 
cavitation  field  ahead  of  the  shock  is  produced  by  interactions  between 
the  liquid  and  container  wall;  the  reaction  required  to  support  the  wave 
is  initiated  in  bubbles  compressed  by  the  shock.  This  model  is  called 
the  cavitation  model  of  LVD. 
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Propagation  of  LVD  is  thus  a  very  complex  phenomenon  involving  the 
formation  of  a  cavitation  field  and  its  subsequent  behavior  under  shock 
compression.  A  detailed  treatment  must  account  for  liquid-container 
wall  interactions  and  also  for  the  dynamics  and  thermodynamics  of  a 
cavitation  field  with  heat  transfer,  mass  transfer,  and  chemical  reaction. 
Here,  however,  it  is  assumed  that  a  detailed  account  of  these  processes 
is  not  required  to  calculate  the  propagation  velocity  of  LVD .  The  pro¬ 
cesses  will  be  modeled  simply  with  well-defined  assumptions  related  to 
their  characteristic  times.  This  approach  is  motivated  by  the  one¬ 
dimensional  C-J  model  for  gaseous  detonation,  which  is  satisfactory  for 
calculating  the  average  velocity  of  self-sustaining  detonation  waves 
without  treating  the  complex  system  of  interacting  transverse  waves 
necessary  for  their  propagation.  As  with  the  C-J  model,  the  validity 
of  assumptions  and  the  model  based  upon  them  must  be  determined  by  the 
comparison  of  calculated  LVD  propagation  velocities  with  those  obtained 
experimentally. 

Hydro thermodynamic  Assumptions 

Our  one-dimensional  treatment  of  LVD  is  based  on  the  following 
assumptions  ordering  the  time  scales  for  cavitation,  mechanical  equilib¬ 
rium,  and  chemical  reaction.  The  cavitation  field  is  assumed  to  be 
formed  quickly  ahead  of  the  main  shock  in  the  liquid,  and  under  shock 
compression  is  assumed  to  attain  mechanical  equilibrium  and  form  a 
steady-state  burning  zone.  The  implications  of  the  assiur.pt ions  arc  as 

follows.  The  rapid  formation  of  the  cavitation  field  necessitates  that 

8 

the  shock  in  the  liquid  be  represented  by  a  Rayleigh  line  passing 
through  the  initial  state  of  the  cavitated  liquid  and  not  the  initial 
state  of  the  virgin  liquid.  The  assumption  about  mechanical  equilibrium 
implies  that  pressure  be  well  defined  at  each  point  of  t>»o  reaction 


zone,  and  specifies  the  initial  state  -for  burning  as  the  initial  shocked 
state  produced  in  the  cavitated  liquid.  Tne  steady-state  assumption 
leads  to  the  conclusion  that  the  reaction  zone  lies  on  the  Rayleigh  line 
representing  the  main  shock  in  the  cavitated  liquid.  An  additional 
assumption  about  the  termination  of  the  reaction  zone  is  required  to 
formulate  a  model  for  calculating  the  propagation  velocity  of  LVD.  We 
will  first  consider"  the  case  when  the  reaction  attains  equilibrium  at 
the  termination  of  the  steady  state  and  then  the  case  when  it  does  not. 

Complete  Equilibrium 

Considerations  based  on  the  properties  of  Hugoniot  curves 

lead  to  the  conclusion  that  the  cavitation  model  for  LVD  gives  a  unique 

propagation  velocity  when  the  reaction  zone  ends  in  a  state  of  complete 

equilibrium.  The  Hugoniot  for  cavitated  liquid  is  shown  schematically 

* 

in  the  pressure-specific  volume  (p-v)  plane  as  OSH  in  Figure  1;  the 

8,9 

detonation  branch  of  the  equilibrium  products  Hugoniot  curve  centered 
at  0  is  shown  as  O^JH'  and  the  Rayleigh  line  tangent  at  its  C-J 
point  J  as  OJS.  The  point  of  intersection  of  the  Rayleigh  line  Oi 
passing  through  the  initial  condition  0  ,  with  slope  -(D/v0)  ,  and 

OSH  represents  the  initial  shocked  state  i  produced  by  a  shock 
propagating  at  velocity  D  in  the  cavitated  liquid  with  initial  volume 
v0 .  Since  the  initial  state  for  burning  is  the  initial  shocked  state, 
the  final  burnt  state  must  lie  on  the  deflagration  branch  of  an  equilibrium 
products  Hugoniot  centered  on  OSH.  We  are  therefore  interested  in  the 
deflagration  branches  of  the  family  of  equilibrium  products  Hugoniot 
curves  centered  on  OSH 

For  convenience,  we  will  consider  the  family  of  deflagration  Hugoniot 

curves  centered  on  OSH  above  S  with  p.  >  p  .  A  member  of  tins  family 

is 

is  shown  schematically  in  Figui’O  1  as  0.J,0.,>  When  p.  >  p  ,  we  will 

ii  is  ; 

*  Figures  are  shown  at  the  end  of  this  appendix.  i 

I 
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show  that  the  points  of  intersection  of  the  Rayleigh  line  Oi  with 

the  detonation  branch  of  the  products  Hugoniot  curve  centered  at  the 

initial  conditdor  0  lie  on  the  deflagration  branch  of  the  products 

Hugoniot  curve  centeved  at  i.  The  points  of  intersection  of  Oi  and 

O'JH7  are  shown  in  Figure  1  as  1  and  2.  It  is  convenient  first  to 

relate  the  states  0,  i,  and  2  by  expressing  the  conditions  that  i 

and  2  lie  on  the  Hugonioi  curves  centered  at  0  and  also  on  the 

Rayleigh  lino  02i.  Let  o  denote  specific  internal  energy,  then 

(o,  p,  v)  states  on  a  Hugoniot  curve  are  related  by  the  Hugoniot  equa- 
8 

tion  expressing  the  balance  of  mass,  momentum,  and  energy  across  a 
shock  discontinuity.  The  Hugoniot  equation  for  a  Hugoniot  curve 
centered  at  0(eo,PoiV0)  is 


2(e  -  e0)  =  (p  +  Po)(v0  -  v)  ( 1) 


and  the 

conditions 

for  i 

to 

be  on  OSH 

and  2 

to  be  on 

o'jh'  follow 

f  rom  Eq , 

.  (1)  as 

2(e.  - 

'  Co) 

=  (Pt  + 

O 

O 

1 

v.) 

i 

(2) 

to 

/-N 

o 

ft) 

1 

■  Go) 

=  ( Pa  + 

Po)(v0  - 

v2) 

(3) 

The  condition  for  these  states  to  bo  on  the  Rayleigh  line  Oi  follows 

8 

from  the  Rankine-Hugoniot  jump  conditions  as 

=  _  (4) 

Vo  -  v2  v0  - 

The  combination  of  Eqs.  (2),  (3),  and  (d)  to  eliminate  e0 ,  p0 , 
and  v0  leads  to  the  equation 

2(0;,  -  e^)  =  (p3  +  p^Cv  -  v2)  (5) 
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relating  the  states  i  and  2,  and  a  similar  argument  loads  to  the 
equation 

2(ex  -  e.)  =  (pi  +  p.)(v.  -  Vj )  (6) 

1  ii 

relating  the  states  i  and  1.  Equations  (5)  and  (6)  show  that  states 
2  and  1  satisfy  the  Hugoniot  equation 

2(e  -  ej  =  (p  +  p  )(v  -  v)  (7) 

for  the  deflagration  branch  of  the  products  Hugoniot  0  J(0  1  centered 

i  i 

at  i.  Thus  when  p.  >  p  ,  the  deflagration  branch  of  the  products 

1  s 

Hugoniot  centered  at  i  and.  the  detonation  branch  of  the  products 
Hugoniot  centered  at  0  intersect  at  the  points  where  the  Rayleigh 
line  connecting  0  and  i  intersects  the  detonation  branch  centered 
at  O. 

It  is  clear  from  Figure  1  that  the  intersection  points  1  and  2 

approach  each  other  as  p  approaches  p  and  coincide  at  J  when 

i  s 

p  =  p  .  Moreover,  when  p  >  p  the  Rayleigh  line  iJ /  tangent  to 
is  is 

0  J  /0  '  intersects  the  p  =  0  axis  at  a  point  to  the  left  of  the 
i  i 

initial  condition  0,  When  p^  <  p  ,  the  deflagration  branch  of  the 
Hugoniot  centered  at  1  does  not  intersect  the  detonation  branch  of 
the  Hugoniot  centered  at  i  ,  and  its  tangent  Rayleigh  line  intersects 
the  p  =  0  axis  at  a  point  lying  to  the  right  of  0.  A  unique  situa¬ 
tion  arises  when  p  =  p  >  for  then  the  Hugoniots  centered  at  0  and 
i  s 

i  are  tangent  to  the  Rayleigh  lino  OS  at  J  ,  and  the  C-J  points  of 
both  curves  coincide.  This  condition  is  the  only  one  compatible  with 
our  assumption  of  a  steady  state  terminating  in  complete  equilibrium 
because  in  this  case  the  Rayleigh  line  must  be  tangent  to  the  deflagra¬ 
tion  branch  of  the  Hugoniot  centered  on  the  initial  shocked  state  and 
must  also  pass  through  the  initial  state  of  the  cavitated  liquid. 
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Thus  conditions  for  a  steady  state  define  a  unique  propagation  velocity 
in  the  simplest  case  when  the  reaction  attains  equilibrium.  The  propa¬ 
gation  velocity  for  LVD  with  the  reaction  attaining  equilibrium  in  the 
cavitation  model  is  determined  by  the  C-J  condition  on  the  detonation 
branch  of  the  equilibrium  Hugoniot  centered  on  the  cavitated  liquid. 

To  test  the  model  with  a  realistic  equation  of  state,  the  TIGER 
code  developed  at  Stanford  Research  Institute  for  the  Ballistic  Research 
Laboratories  under  Contract  No.  DA-04-200-AMC-3226(X) ,  was  used  to  calcu¬ 
late  the  detonation  parameters  of  ethyl  nitrate  for  different  degrees 
of  cavitation.  The  internal  energy  of  the  cavitated  liquid  ahead  of  the 
shock  was  assumed  to  be  the  same  as  that  of  the  uncavitated  liquid.  The 
results  of  these  calculations  are  shown  in  Table  1.  Calculations  with 
the  BKW  and  the  virial  equation  of  state  do  not  agree  with  experiment, 
even  when  the  volume  of  the  cavitated  liquid  greatly  exceeds  that  expected 
experimentally.  Propagation  velocities  calculated  with  the  ideal  gas 
equation  of  state  agree  with  experiment,  but  the  detonation  pressures 
lie  in  the  kilobar  region  where  this  equation  of  state  is  inapplicable. 

The  results  of  these  calculations  lead  to  the  conclusion  that  the 
assumption  of  complete  reaction  is  not  valid.  Attention  was  therefore 
given  to  a  study  of  the  cavitation  model  when  the  reaction  does  not 
attain  equilibrium. 

Incomplete  Equilibrium 

We  will  use  a  reaction  coordinate  X  to  denote  the  fraction  of 
unburnt  liquid  when  the  reaction  does  not  attain  equilibrium.  Then 
X  =  0  in  the  former  case  when  the  reaction  attains  equilibrium  and  all 
the  liquid  is  consumed,  and  X  =  1  when  there  is  no  reaction.  Our 
treatment  of  partial  equilibrium  is  based  on  the  assumption  that  the 
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reaction  is  frozen  with  the  products  in  chemical  equilibrium.  In  this 
case  a  value  of  X  defines  an  equilibrium  products  Hugoniol  curve  for 
incomplete  reaction — a  frozen  Hugoniot — centered  on  the  initial  condition 
of  the  cavitated  liquid.  And  the  values  of  X  in  the  range  0  <  X  <  1 
define  a  family  of  frozen  Hugoniot  curves  lying  between  the  Hugoniot 
curve  of  the  cavitated  liquid  and  the  equilibrium  products  Hugoniot  for 
complete  reaction. 

It  follows  from  the  previous  treatment  of  complete  equilibrium  that 
LVD  will  be  represented  by  a  Rayleigh  line  tangent  to  the  detonation 
branch  of  a  frozen  Hugoniot  centered  on  the  initial  state  o  of  the 
cavitated  liquid.  An  equation  of  state  for  nitromethane  was  incorporated 
into  the  TIGER  code  so  that  detonation  parameters  could  be  calculated 
with  a  realistic  thermodynamic  description  of  the  explosive  and  the 
products  of  reaction.  Detonation  parameters  wore  calculated  as  a  function 
of  X  to  determine  if  LVD  could  be  modeled  with  the  cavitation  model, 

10 

A  complete  equation  of  state  of  the  liquid  explosive  must  be 
incorporated  into  the  TIGER  code  to  compute  detonation  parameters  when 
the  reaction  does  not  proceed  to  equilibrium.  A  knowledge  of  the  incom¬ 
plete  pressure-volume-temperature  (p-v-T)  equation  of  state  and  the  varia¬ 
tion  of  specific  heat  at  constant  pressure  C  along  the  atmospheric 
isobar  is  rcuired  to  introduce  a  complete  equation  of  state  into  the 
TIGER  code.  The  thermodynamic  functions  required  to  compute  a  thermo¬ 
dynamic  state  are  calculated  by  integrating  thermodynamic  identities 
along  isotherms  from  the  atmospheric  isobar.  The  (p-v-T)  relationship 
used  in  the  present  work  was  based  on  the  equation  of  state  for  liquids 
formulated  previously  for  shock  temperature  calculations.  In  this 
description  of  liquids,  Op/9T)  is  assumed  to  be  constant  and  the 
specific  heat  at  constant  volume  Is  assumed  to  be  a  function  of 

temperature,  i.e.  C  (T) .  Nitromethane  was  chosen  as  a  typical  liquid 
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explosive  because  C  (T)  and  the  shock  temperature  along  the  Hugoniot 

V  11 

curve  were  known  from  the  previous  work.  The  standard  Hugoniot  centered 

o 

on  the  initial  state  (p  =1  atmosphere,  T0  =  298  K,  and  v0  =  0.884 

o 

cc/g)  in  the  (p-v)  plane  was  constructed  with  the  universal  Hugoniot 

for  liquids  U  =  aiC0  +  a2u  where  aj.  and  a2  are  constants;  and  U,  u, 

and  c0  denote  shock  velocity,  particle  velocity,  and  sound  speed  in 

* 

the  initial  condition. 

Let  b  =  Op/ST)  ,  a  =  l/v0Ov/5T)  ,  and  let  the  subscript 

v  p=l 

h  denote  states  on  the  standard  Hugoniot  curve  shown  schematically  as 
03h  in  Figure  2.  The  (p-v-T)  equation  of  state  in  the  region  v  v0 
spanned  by  the  Hugoniot  curve  can  be  written  as 

p  =  p  (v)  +  b[T  -  T(v)],  v  :£  v0  (8) 

h  h 

where 

Pj(v)  =  p0  +  (axc0)2  (v0  -  vh)/[v0  -  a2(v0  -  v^]2  (9) 


and  T^(v)  is  calculated  by  integrating  the  differential  equation 
governing  temperature  along  the  Hugoniot  curve.  The  (p-v-T)  equation 
of  state  was  extended  to  the  right  of  v0  by  assuming  that  Ov/ST) 
was  constant  and  therefore  equal  to  0Cvo  along  the  atmospheric  isobar 
p  =  p0  =  1  when  v  S  v0 .  With  this  assumption  the  (p-v-T)  equation 
of  state  in  the  region  v  s  v0  can  be  written  as 

p  -  p0  =  b[(T  -  T0)  -  l/O(v/v0  -  1)],  v  2:  v0  .  (10) 


Additional  work  on  the  universal  Hugoniot.  curve  for  liquids  is 
presented  in  this  report  in  Appendix  C. 


/ 

/ 


) 
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Equations  (8)  and  (10)  have  continuous  first  derivatives  where 
they  are  patched  together  along  the  v0  isochore  042.  In  performing 
an  integration  along  an  isotherm,  say  543  (Fig,  2),  Eq.  (10)  is  used 
until  v  =  v0  at  4,  and  then  Eq.  (8)  is  used  along  43  for  v  5  v0 . 

Equations  for  the  specific  Gibbs  free  energy  G,  the  specific 
entropy  S  ,  the  specific  enthalpy  11,  and  C  were  derived  with  Eqs. 

p 

(8)  and  (10)  and  incorporated  into  the  TIGER  code.  But  only  the  equa¬ 
tions  for  G  will  be  given  here  because  the  equations  for  the  other 
state  variables  can  be  calculated  from  them  with  thermodynamic  identities. 
For  v  s  v0  ,  the  equation  for  G  is 

G(T , v)  =  G(T)  -  b/2  av0[v2  -  (v(T))R]  (11) 

with 

T  T 

G(T)  =  G0  -  So(T  -  T0)  -  T  [  (C  /T)  liT  +  I"  C  (IT  (12) 

JT0  P  To  P 

and 

v(T)  =  v0(l  +  a(T  -  T0))  .  (13) 

For  v  £  v0  ,  it  is 

G(T,v)  =  G(T,v0)  +  pv  -  p(T,v0)v0-  bT(v  -  v0)  -  I  +  bIT  (14) 

where  G(T,v0)  and  p(T,v0)  are  obtained  by  setting  v  =  v0  in 

Eqs,  (8)  and  (11),  and  I  =  p  dv  and  I  =  T  dv  denote  the 

p  J  h  T  J  h 

Vo  v0 

integrals  of  shock  pressure  and  temperature  along  the  Hugoniot  curve. 

The  SPECIAL  routine,  developed  at  Stanford  Research  Institute  for 
Picatinny  Arsenal  under  Contract  DAAA21-71-C-0454  to  treat  nonideal 
detonation,  was  used  to  calculate  detonation  parameters  of  cavitated 
nitromethane  as  a  function  of  the  unreacted  fraction  X  with  different 
equations  of  state  for  the  reaction  products.  Most  of  the  calculations 


were  performed  with  the  virial  equation  of  state  because  we  are  inter¬ 
ested  in  detonation  pressures  in  the  region  of  10  kilobars.  The  virial 
equation  of  state  used  in  the  TIGER  code  contains  the  second  virial 
coefficient  B  and  the  third  virial  coefficient  C.  The  results  of 
the  calculations  performed  without  the  third  virial  coefficient  (C  =  0) 
for  a  10  percent  increase  in  initial  volume  are  given  in  Table  2,  and 
the  results  of  the  calculations  performed  with  the  BKIV  and  ideal  gas 
equations  of  state,  and  the  virial  equation  with  C  ?■  0  aro  also  given. 
Because  of  convergence  problems  in  the  SPECIAL  routine,  the  calculations 
for  a  10  percent  increase  in  initial  volume  were  limited  to  values  of 
X  in  the  range  0  <  X 0.625,  and  the  case  of  a  20  percent  increase  in 
the  initial  volume  was  considered. 

The  results  of  the  calculations  performed  for  a  20  percent  increase 
in  initial  volume  using  the  virial  equation  with  C  =  0  are  given  in 
Table  3.  The  virial  equation  with  C  =  0  is  used  for  all  values  of  X 
for  convenience  because  the  SPECIAL  routine  was  formulated  to  calculate 
detonation  parameters  along  the  frozen  C-J  locus  from  the  equilibrium, 

X  =  0,  C-J  state.  The  calculated  detonation  pressures  and  velocities 
for  X  in  the  range  0  <  X  S  0,75  are  plotted  as  a  function  of  X  in 
Figure  3.  For  values  of  X  >  0.75  the  graphs  were  extrapolated  smoothly 
to  satisfy  the  boundary  conditions  that  the  pressure  approaches  one 
atmosphere  and  the  velocity  approaches  the  sound  speed  in  the  cavitated 
liquid  as  X  approaches  one.  The  sound  speed  in  cavitated  liquid  was 
estimated  to  be  1,88  mm/(isoc  when  the  increase  in  initial  volume  was  10 
percent  and  to  be  1.08  mm/p,sec  when  the  increase  was  20  percent. 

It  follows  from  Figure  3  that  the  calculations  are  not  completely 
definitive  because  experimental  values  of  LVD  velocity  and  pressure 
correspond  to  values  of  X  w  0.8  and  lie  in  the  region  where  the  results 
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of  the  calculations  had  to  be  extrapolated.  Further  examination  of 
Figure  3,  however,  shows  that  the  curves  for  a  10  and  20  percent  increase 
in  initial  volume  lie  very  close  together  in  the  10-kilobar  region  and 
leads  to  the  conclusion  that  a  surprisingly  small  amount  of  reaction  is 
required  to  support  LVD.  Important  features  of  the  cavitation  model 
with  incomplete  reaction  can  therefore  be  stated  as  follows: 

(1)  Only  reaction  of  about  20  percent  of  the  liquid  is 
required  to  support  LVD. 

(2)  LVD  parameters  do  not  depend  strongly  on  the  initial 
degree  of  cavitation. 

Conclusions 

It  is  clear  that  additional  work  is  required  to  formulate  a  satis¬ 
factory  method  for  calculating  the  propagation  velocity  of  LVD.  This 
work  should  be  concerned  mainly  with  the  degree  of  incomplete  reaction. 

An  attempt  should  be  made  to  establish  a  condition  for  the  termination 
of  the  reaction  zone  so  that  there  is  a  unique  propagation  velocity  for 
the  cavitation  model  with  incomplete  reaction.  A  consideration  of  bubble 
dynamics  with  heat  and  mass  transfer  and  chemical  reaction  should  be 
undertaken  to  try  and  determine  a  criti*”  lity  condition  for  termination 
of  reaction  and  a  means  of  calculating  a  value  of  X  at  the  end  of  the 
reaction  zone.  LVD  velocities  for  different  explosives  calculated  with 
values  of  X  could  then  he  used  to  test  the  validity  of  the  cavitation 
model.  The  detonation  products  of  a  liquid  explosive  undergoing  first 
I1VD  and  then  LVD  should  be  collected  and  compared  to  test  if  the  amount 
of  reaction  predicted  by  the  present  work  with  the  cavitation  model  for 
LVD  is  reasonable. 
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*  Detonation  in  virgin  material  (no  cavitation) 
**  Third  virial  coefficient  set  equal  to  zero. 


Table  2 

CALCULATED  DETONATION  PARAMETERS  FOR  NITROMETHANE 
WITH  THE  INITIAL  VOLUME  INCREASED  TEN  PERCENT 


Equa  t ion 
of 

State 


Virial 
(C  =  0) 


Virial 
(B  and  C) 


Ideal  Gas 


Unreacted 

Fraction 

(X) 

Detonation 

Velocity 

(mm/(j.sec) 

C-J 

Pressure 
(atm.  x  103) 

C-J 

Volume 
(  cc/g) 

Adiabatic 

Index 

(x> 

0.0 

4.04 

54.9 

0.650 

2.02 

1/16 

3.91 

53.9 

0.635 

1.88 

1/8 

3.79 

51.4 

0.629 

1.83 

3/16 

3.66 

48.6 

0.625 

1.80 

1/4 

3.53 

45.6 

0.623 

1.78 

3/8 

3.29 

39.5 

0.622 

1.78 

1/2 

3.05 

33.7 

0.625 

1.80 

5/8 

2.80 

27.5 

0.635 

1.88 

0.0 

6.49 

118.7 

0.702 

2.59 

1/4 

5.55 

81.9 

0.717 

2.82 

1/2 

4.47 

50.8 

0.729 

2.99 

0.0 

5.20 

79.5 

0.691 

2.45 

1/8 

4.70 

65.1 

0.689 

•J.44 

1/4 

4.24 

53.7 

0.685 

2.39 

3/8 

3.82 

44.5 

0.679 

2.32 

0.0 

2.32 

24.8 

0.533 

1.22 

1/4 

2.62 

33.6 

0.503 

1.07 

1/2 

2,44 

28.4 

0.515 

1.13 
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The  universal  Hugoniot  curve1-3  used  previously  to  calculate  shock 
temperature  in  liquid  explosives4'5  is  not  valid  below  20  kilobars  be¬ 
cause  it  was  constructed  from  shock  wave  data  obtained  at  higher  pres¬ 
sures.  Moreover,  it  does  not  satisfy  conditions  defined  in  the  initial 
state  by  the  Rankine-Hugoniot  jump  .conditions.  The  object  of  the  pres¬ 
ent  work  was  to  extend  the  Hugoniot  curve  to  the  initial  state  and  thus 
improve  the  thermodynamic  description  of  liquids  in  the  1-bar  to  20- 
kilobar  region.  The  objective  was  attained  by  modifying  the  form  of  the 
Hugoniot  curve  to  satisfy  the  initial  conditions,0  calibrating  it  with 
static  pressure  data  for  water,7'8  and  checking  its  validity  with  dynamic 
shock  wave  data  for  glycerine9  and  carbon  tetrachloride. 10 

The  Hugoniot  curve  defines  the  locus  of  shocked  states  obtainable 
from  a  given  initial  condition.  It  is  obtained  experimentally11  from  a 
determination  of  the  states  produced  by  constant  velocity  shocks  propa¬ 
gating  at  different  velocities.  Experimental  determination  of  a  Hugonoit 
is  usually  very  expensive  and  time-consuming,  and  requires  considerable 
amounts  of  materials  which  are  destroyed  in  the  process.  There  is  there¬ 
fore  a  need  for  estimating  Hugoniots  from  easily  measured  physical 
properties. 

Following  earlier  work  by  Gibson  and  coworkers,1  Wool-folk  and  Amster  2 
and  Voskoboinikov,  Afanasenkov,  and  Bogomolov  3  have  shown  that  the  Hugo- 
noits  of  liquids  could  be  represented  by  a  single  normalized  plot  of  the 
form 


U/co  =  a,  +  (a2u/co) 


(1) 


whore  U  is  the  shock  velocity,  u  is  the  particle  velocity,  c  is  the  sound 
velocity,  a,  and  a2  are  constants,  and  subscript  zero  denotes  the  initial 
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4 


1.62 


state  at  a  pressure  p  of  1  bar.  Eqn.  (1)  with  a,  =  1.37  and  a2  = 

fitted2  the  experimental  data  which  were  in  the  range  u/c  =  0.5  to  2.5, 

o 

corresponding  to  shock  pressures  in  the  range  of  20  to  150  kbar. 

The  main  problem  with  Eqn.  (1)  is  that  it  does  not  satisfy  the  bound¬ 
ary  condition  U  =  c  at  u  =  0.  Eqn.  (1)  therefore  cannot  be  used  in  the 
o 

region  from  1  bar  to  20  kbar.  This  low-pressure  region  is  of  interest  be¬ 
cause  of  its  importance  in  low-velocity  detonations. 

Jacobs0  has  suggested  that  an  additional  term  (1-aj)  exp[-a3u/co]  with 
a3  constant,  be  added  to  Eqn.  (1)  so  that  the  boundary  condition  may  be  met. 
The  form  of  the  "universal"  Hugoniot  would  then  be 

U/c  =  1.37  -  (0.37)  exp[-a,u/c  ]  +  1,62  u/c  (2) 

which  reduces  to  U/cq  =  1  at  u  =  0.  In  Fig.  1  the  experimental  data  are 
shown  along  with  three  calculated  curves  that  correspond  to  a3  =  «,,  5,  and 
1.  Figure  1  indicates  that  the  value  of  a3  should  be  <.  5.  An  express  !on 
for  a3  was  derived  from  thermodynamic  identities  relating  the  Hugoniot 
curve  and  the  isentrope.  The  constant  a3  was  evaluated  using  echo-sounding 
data  for  water. 

The  identity  for  the  initial  slope  of  tiie  Hugoniot  in  the  (U-u)  plane 


o 


where  v  denotes  specific  volume  and  s  denotes  specific  entropy,  and  the 
thermodynamic  identity  for  an  isentrope 
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(4) 


2c 


were  used  to  determine  the  constant  a3  in  Eqn,  (2).  Differentiating  Eqn. 
(2)  with  respect  to  u  gives  the  slope  of  the  Hugoniot  as 

dU/du  =  1.62  +  0.37  a3  expC-a3u/cQ)  (51 

and  the  relationship  between  the  initial  slope  of  the  Hugoniot  and  a3  is 
obtained  as 


0.37  a_  =  ((dU/du)  ,  -  1.62) 
3  '  p=l 


(6) 


by  setting  u  =  0  in  Eqn.  (5).  The  calculation  of  (52 p/dv2l  at  p  =  1  bar 

s 

with  sound  velocity  data  is  then  sufficient  to  evaluate  a3.  The  identity 


(7) 


was  used  in  the  evaluation  of  a3  because  it  is  convenient  to  estimate  the 

derivatives  (Sc/5v)  and  (Qv/Sp)  with  experimental  data.  Data  from  echo- 
p  c 

sounding  tables  for  pure  water,7  together  with  p-v-T  data,8  were  used  to 

construct  a  graph  of  the  variation  of  sound  speed  with  specific  volume  at 

pressures  of  1,  50,  100,  and  150  bar  (see  Figure  2).  At  p  =  1  bar,  T  = 

25° C,  and  c  =  1493  m/sec,  we  estimate  (dc/dv)  .  from  the  graph  as 
o  '  '  p= 1 


(  ~)  =  1.027  x  10°  g/(scc  cm2) 

VAV  i=l 


c  -:j 


If 


Similarly,  from  the  graph,  at  constant  c  =  1493  m/sec,  the  following  val¬ 
ues  of  p  and  v  were  estimated: 


p 

V 

■Lv 

(bar) 

(cc/g) 

(cc/g) 

1 

1.00295 

0.00325 

50 

0,99970 

0.00290 

100 

0.99680 

0.00275 

150 

0.99405 

Extrapolating  to  p  =  1,  we  obtain 


(I?)  =(tp)  = '7  »  10"”  e) 


Substitution  of  c  =  1.493  X  105  cm/sec,  v  =  1.00295  cm3/g,  (dv/dp)  = 

o 

-7  x  10-11  cm5/(dyne  g),  and  Oc/Sv)^^  =  1.027  x  10G  g(sec  cm2)  into 

Eqn.  (7)  gives  (dc/dv)  =  5.66  x  105  g/(sec  cm2).  Evaluating  (d2p/dv2)  = 

s  s 

2.12  X  1011  sec2/(g3  cm"7)  with  Eqn.  (4)  and  (dU/du)^^  =  2.40  with 
Eqn.  (3)  gives  a3  =  2.1  by  substitution  in  Eqn.  (5). 

The  region  of  Fig.  1  close  to  the  origin  has  been  replotted  as  Fig.  3, 
using  Eqn.  (2)  and  values  of  a3  =  ra,  10,  5,  2,  1,  and  0.1.  Also  plotted 
in  Fig.  3  are  experimental  Hugoniot  points  for  glycerin  obtained  by  Erlich9 
and  for  carbon  tetrachloride  obtained  by  Lysne. 10  Until  further  data  are 
available,  we  conclude  that  a  value  of  a3  =  2  is  consistent  with  results 
of  the  calculations  on  water  and  with  the  experimental  data  on  glycerin 
and  carbon  tetrachloride.  In  other  words,  Eqn.  (2)  becomes: 


U/cQ  =  1.37  -  0.37  exp[-2u/cQ]  +  1.62  u/c^ 


(8) 


It  is  concluded  that  the  modified  universal  Hugoniot  curve  formulat¬ 
ed  in  this  note  improves  the  thermodynamic  description  of  liquids  in  the 
1-bar  to  20-kbar  region. 
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Introduction 


The  hot  spots  caused  by  adiabatic  heating  during  the  collapse  of 

1 

bubbles  or  cavities  wore  proposed  by  Bowden  and  coworkers  in  1947  ns 

a  mechanism  for  the  initiation  of  explosion  in  highly  energetic  liquids. 

This  initiation  mechanism  was  proposed  to  explain  the  high  sensitivity 

of  liquid  or  gelatinous  explosives  to  mechanical  impact.  As  a  result  of 

extensive  subsequent  investigations,  the  effect  of  hot  spots  on  the  shock 

sensitivity  of  liquid  explosives  is  well  understood.  The  early  experi- 
2 

monts  by  Yoffo  sipported  the  view  that  initiation  is  duo  to  compressive 
heating  of  entrapped  gas  or  vapor  spaces. 

The  detailed  mechanism  of  hot-spot  initiation  by  bubbles  in  some 

cases  may  be  due  to  processes  other  than  heating  by  adiabatic  compression. 

3 

For  example,  studies  by  Campbell  ot  al.  show  that  initiation  may  take 

place  by  shock  interactions  caused  by  both  solid  inhomogeneities  and 
4 

bubbles.  Mader,  in  a  two-dimensional  numerical  computation  of  the 

reactive  flow  around  cylindrical  and  rectangular  cross  section  rods  or 

cavities  placed  parallel  to  the  shock  front,  has  given  added  support  to 

the  conclusions  reached  in  those  experiments.  Heating  sufficient  to 

cause  initiation  by  cylindrical ly  symmetric  shock  interaction  with  a 

5 

spherical  bubble  was  also  computed  by  Evans  ot  al.  Hot-spot  ignition 

caused  by  shock  interaction  or  focusing  at  a  bubbl  is  of  course 

increasingly  effective  with  higher  shock  strengths.  On  the  other  hand, 

it  is  also  conceivable  that  the  Munroc  jot  effect  could  be  responsible 

for  bubble  and  cavity  initiation  at  higher  shock  strengths.  This 

6 

mechanism  was  proposed  by  Bowden  and  subsequently  microjet  ignition 

7 

was  observed  experimentally  by  Watson  and  Gibson  and  under  different 

8 

conditions  by  Bowden  and  McOnlo.  At  lower  shock  strengths,  however, 

9 

Gibson  ot  al.  have  shown  that  the  velocity  of  the  jet  is  not  sufficient 
to  cause  shock  initiation  as  the  jet  collides  with  the  opposite  side 


of  the  bubble. 
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Thus,  hot-spot  initiation  by  a  bubble  clue  to  both  shock  focusing  and 
microjetting  are  generally  of  interest  l'or  shock  strengths  in  excess 
of  about  50  kbar. 


The  important  role  of  bubble  hot-spot  ignition  in  the  initiation 
and  propagation  of  low-velocity  detonation  (LVD),  however,  has  recently 
stimulated  renewed  interest  in  hot-spot  formation  at  lower  shock  over¬ 
pressures  in  the  1  to  10  kbar  range.  The  .most  plausible  model  of  LVD 

10 

is  a  cavitation  model  recently  suggested  by  Watson  et  al.  and  inde¬ 
pendently  by  Voskoboinikov  et  al.^  In  the  cavitation  model,  precursor 
container  wall  waves  cavitate  the  liquid  ahead  of  a  liquid  shock;  the 
cavitation  bubbles  are  subsequently  compressed  by  the  liquid  shock 
causing  hot-spot-initiated  reaction;  and  the  resulting  energy  release 

in  turn  supports  the  liquid  shock.  This  model  has  been  further  verified 
12 

by  Watson  who  has  photographed  the  cavitation  ahead  of  the  reaction 

13 

zone  and  by  Watson  et  al.  who  more  recently  have  reported  tentative 

measurements  of  precursor  wall  waves  capable  of  producing  cavitation 

ahead  of  the  reaction  zone.  In  addition,  the  existence  of  individual 

burning  sites  rather  than  a  continuous  reaction  front  has  been  observed 
12 

by  Watson  in  x-ray  photographs  of  the  LVD  reaction  wave,  and  Gibson 
14 

et  al.  have  experimentally  verified  the  dependence  of  bubble  ignition 
on  the  strong tli  of  the  initiating  shock. 

'Hie  cavitation  model  also  explains  the  relatively  long  dark  region 

12 

observed  between  the  shock  front  and  the  luminescent  reaction  zone. 

Thus,  the  cavitation  model,  though  it  involves  an  incompletely  under¬ 
stood  interaction  between  the  container  and  the  liquid  explosive,  the 
formation  of  cavitation  and  subsequent  hot-spot  reaction,  and  as  yet 
unknown  criteria  for  stability,  appears  to  offer  the  most  plausible 
mechanism  for  the  propagation  of  LVD. 
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c 


Background 


The  dynamics  of  vapor  bubbles  has  been  a  subject  of  extensive 

theoretical  and  experimental  investigation  for  a  number  of  years.  In 

15-17 

a  sequence  of  early  papers,  Plesset  and  Zwick  calculated  the  rate 

of  growth  of  vapor  bubbles  in  slightly  superheated  steam  under  constant 
external  pressure.  Their  calculations  assume  that  the  temperature 
variations  in  the  liquid  are  appreciable  only  in  a  thin  thermal  boundary 
layer  near  the  bubble  wall.  The  predicted  bubble  radius  R(t)  is 

t 

o 

asymptotically  proportional  to  t“  when  R  is  large  enough  so  that  surface 

tension  is  negligible,  a  result  that  was  confirmed  experimentally  by 
18 

Dergarabcdian.  More  recent  work  has  extended  the  work  of  Plesset 
15-17  19-21 

and  Zwick  to  more  general  cases.  Most  of  the  work  on  the 

dynamics  of  vapor  bubbles  has  been  aimed,  however,  at  boiling  and 
related  phenomena.  Of  more  interest  hero  is  the  collapse  of  a  bubble 
leading  to  high  interior  temperatures  causing  local  initiation  of  re¬ 
action  which  may  ultimately  ignite  the  bulk  material  leading  to 
explosion. 

Theoretical  studies  of  hot-spot  initiation  that  are  related  to 

22 

the  adiabatic  compression  of  a  bubble  have  been  carried  out  by  Zinn, 

25  2d 

Ivnxg,  and  Gill.  Those  are  computations  of  ignition  delay  time 
which  arc  applicable  to  a  bublJlo  with  a  stationary  radius  and  with  no 
mass  transfer  at  the  bubble  wall.  In  the  present  computations,  however, 
wo  are  interested  in  the  temperature  rise  and  the  onset  of  reaction 
or  ignition  in  a  collapsing  vapor  bubble.  We  include  here  the  effects 
of  both  heat  and  mass  transfer  (vaporization)  at  the  bubble  wall,  and 
formulate  the  problem  in  Lagrangian  coordinates  that  move  with  the 
bubble  wall  during  collapse.  It  is  assumed  that  both  the  vapor  and 
liquid  are  aonviscous  and,  prior  to  collapse,  they  have  the  same 
temperature  T  .  The  interior  of  the  vapor  bubble  is  assumed  to  be 
uniform,  an  assumption  that  is  justified  provided  that  the  thermal 
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diffusivity  of  t!io  vapor  is  much  greater  than  that  of  the  liquid  and 

that  the  speed  of  the  bubble  wall  motion  is  much  loss  than  the  sound 

speed  in  the  vapor.  The  first  condition  is  certainly  true  for  most 

substances  for,  although  the  thermal  conductivity  of  the  liquid  is 

about  ten  times  greater  than  that  of  the  vapor,  the  thermal  diffusivity 

of  the  vapor  is  still  much  greater  than  that  of  the  liquid.  The  second 

25 

condition  has  been  discussed  by  Hickling  in  connection  with  the 

collapse  of  a  bubble  of  inert  gas  in  a  liquid.  Hickling  found  that  the 

sound  speed  in  the  gas  was  indeed  much  higher  than  the  velocity  of  the 

bubble  wall  for  initial  bubble  radii  of  R  =  10  ^  -  10  cm.  and 

o 

overpressures  of  about  4  atm.  Under  such  conditions,  the  pressure  in 
the  vapor  bubble  remains  approximately  uniform  because  the  effects  of 
disturbances  at  the  wall  will  have  time  to  bo  propagated  throughout 
the  bubble  interior.  It  is  clear  that  during  the  early  stages  of 
collapse  the  assumption  of  uniform  bubble  interior  will  certainly  bo 
valid  in  the  present  calculation.  Further,  to  keep  the  problem  tractable, 
a  spherical  bubble  will  be  assumed  oven  though  shock  compression  in 
general  produces  asymmetric  collapse  and,  in  addition,  it  is  known  that 
spherical  bubbles  arc  generally  unstable  during  the  latter  stages  of 
collapse. 
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Formulation  of  the  Problem 


A  very  general  formulation  of  the  equations  governing  the  motion 

26 

of  a  bubble  in  a  liquid  was  given  by  Hsieh,  as  well  as  simplified 
equations  for  a  spherical  bubble  with  a  uniform  interior.  With  the 


above  assumptions  and  the  additional  assumptions  that: 


(1)  No  body  force  is  present, 

(2)  There  is  no  translational  motion  of  the  bubble  relative 
to  the  liquid; 

(3)  Viscosity  coefficients  are  zero  throughout; 

(4)  The  liquid  is  incompressible; 

(5)  The  equation  for  the  temperature  in  the  liquid 
T  (r,  t)  can  bo  written  as 

Jj 


*1 

or 

\r  _ 

3T  + 

1L  Or 

l 


(1) 


where  r  is  the  radial  coordinate  from  the  bubble  center  and  ft.  is 

a 

the  thermal  diffusivity  of  the  liquid,  which  will  be  assumed  constant. 
Since  the  liquid  is  treated  as  incompressible,  and  the  density  of  the 
vapor  0v  is  generally  much  loss  than  that  of  the  liquid,  from  the 
continuity  requirement  the  liquid  velocity  may  bo  directly  related  to 
the  bubble  radius  K(t)  as  follows: 


V 


i 


K“R 

2 


r 


(2) 


That  is,  the  contribution  of  mass  transfer  (vaporization)  at  the  bubble 
surface  to  R  is  neglected.  Note  that  the  assumption  of  uniform  bubble 
interior  requires  that  the  quantities  inside  the  bubble  p  ,  T  ,  P  arc 

V  V  V 

functions  of  time  only.  The  initial  condition  will  be  that  at  t  =  0 

a  bubble  of  radius  It  is  at  rest,  with  die  vapor  inside  at  initial 

o 
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density  and  pressure  of  o  ,  P  ,  respectively.  The  assumption  of 

vo  vo 

uniform  initial  temperature  T  for  the  system  requires  that 

o 

T  (0)  =  T  (r,  0)  =  T  .  (3) 

V  J6  O 

Further,  the  boundary  condition  at  the  bubble  wall  requires  that 

T  (t)  »  T  .(0, t)  (4) 

v  l 

and 


when  I,  is  the  latent  heat  per  unit  mass,  C  is  the  specific  heat 

v 

at  constant  volume  of  the  vapor,  and  Ag  is  the  thermal  conductivity 

of  the  liquid.  The  quantities  L,  and  C  are  assumed  constant, 

l  v 

as  are  Q,  Z,  and  E  which  denote,  respectively,  the  specific  heat  of 
reaction,  the  frequency  factor,  and  the  activation  energy.  'Hie  last 
term  in  Eq .  (5)  is  the  heat  release  in  the  bubble  duo  to  gas  phase 
reaction  which  is  generally  negligible  until  T  roaches  the  order  of 
E/R. 

Since  the  (uniform)  temperature  in  the  vapor  bubble  is  determined 
by  the  temperature  of  the  liquid  at  the  bubble  wall,  it  follows  that 
the  (uniform)  pressure  in  the  bubble  is  given  by  the  vapor  pressure 

of  the  liquid  at  the  wall  or  bubble  temperature.  We  will  therefore 
make  the  additional  assumption  that  the  vapor  pressure  is  the  bubble 
wall  and  that  the  bubble  contains  pure  vapor.  Thus,  the  possibility  that 
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absorbed  gas  in  the  liquid  may  contribute  to  the  interior  pressure  of 
the  bubble  i's  excluded,  and  the  results  will  be  applicable  to  a  pure 
liquid  or  to  cases  in  which  the  partial  pressure  of  absorbed  gases  is 
negligible  in  comparison  with  the  vapor  pressure.  The  equilibrium  vapor 
pressure  requirement  is  certainly  valid  in  the  early  stages  of  collapse 
until  the  rate  of  change  of  the  pressure  in  the  bubble  becomes  comparable 
to  the  kinetic  rate  of  vaporization.  To  enable  us  to  treat  the  problem 
analytically,  we  will  relate  the  vapor  pressure  of  the  liquid,  and  hence 
the  pressure  in  the  bubble,  to  the  temperature  at  the  bubble  surface  by 
use  of  the  Clasuis-Clapeyron  equation 


P  exp 
o 


L 

ft  T 

o 

( 

1 

(6) 


whore  ft  is  the  ideal  gas  constant  and  P  is  the  vapor  pressure  at 

o 

the  initial  temperature  T  .  Further,  the  thermodynamic  quantities  in 
the  vapor  phase  will  be  related  by  the  ideal  gas  equation  of  state 


P  =  ft  o  T 

V  V  V 


(7) 


The  additional  equation  required  to  complete  the  formulation  of  the 
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problem  is  the  differential  equation  for  the  motion  of  the  bubble  wall 

t 

,  9  (Pc=-Pv) 

R  R  +7  (R)2  =  -  - \  t  >  0  ,  (8) 


whore  P^  ,  assumed  constant,  is  the  initial  pressure  in  the  bulk  liquid 
and  corresponds  to  the  pressure  at  largo  distance  from  the  bubble  during 
collapse.  The  initial  conditions  for  Kq .  (8)  are  taken  as 


R(0)  =  R 


(9) 


and 


R (0)  =  0  . 
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Equations  (1),  (2),  (6),  (7),  and  (8)  now  form  a  set  for  the 
dependent  variables  V  (r,t),  T  (r, t),  P  (t),  T  (t),  and  o  (t)  , 

&  X>  V  V  V 

subject  to  the  conditions  of  (3),  (4),  (5),  and  (9). 


Introduce  the  dimensionless  variables 


o 


t  =  t/T 

c 

R  =  R  ( t  )'/R 

o 


T  -  T 


'  l  T 


o 


(10) 


and  the  parameters 


V 

T 

c 


L/RT 

o 

2  i 
(p,r  /p y 

Z  o  0° 


T 

h 


r  2 /a 

o 


l 


T  =  Rt  /QZ  , 
r  o 


(11) 


whore  T  is  the  characteristic  time  for  bubble  collapse  T  is  the 
0  li 

characteristic  time  for  heat  transfer  to  the  bubble  wall,  and  T  is 

r 

the  dimensionless  frequency  factor  for  the  reaction  term  in  Eq.  (5). 


Equation  (1)  can  now  be  written  as 


(12) 


and  the  conditions  (3),  (4),  and  (5)  become  respectively j 


0  (0) 

V 

o 

y'l 

X 

o 

ii 

, - V 

o 

Y, 

CD 

II 

(13) 

S  <t) 

V 

=  0(0,  t)  ,  uo 

(14) 

whore  Y  is  the  ratio  at  specific  heats  of  the  vapor  phase  and  p 

vo 

is  the  initial  density  in  the  vapor  phase. 


The  bubble  equation 

(8) 

can 

be 

written  as 

/  2A\ 

—  / 

,2  P  -  I> 

R  /  Ui\ 

2  / 

dR 

\  ®  V 

(16) 

,  2  / 

+ 

2 

A 

)  P 

\  dt  ' 

\ 

dt 

/  “ 

with  the  initial  conditions 


and 


R(0)  =  l 


dt 


(17) 

(18) 
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Small-Time  Solution 


For  small  time  defined  by  t  «  T  ,  the  bubble  radius  has  the 

c 


form 


:(o) 


*2  . 

1  +  t  /2 


(19) 


so 


a(o) 

dir 


dt 


=  t 


(20) 


Since  9  (0)  =  0,  the  small-time  expansion  of  9  (t)  is  «  1  . 
v  v 

A  similar  statement  may  be  about  about  9^(x,t). 

In  general,  for  small  time,  since  T  «T  ■ for  situations  of 

he 

interest  here,  the  temperature  variation  in  the  liquid  occurs  in  a 
thin  layer  near  the  surface  of  the  bubble  such  that  in  Eq.  (12)  wo  may 

a3  ^ 

neglect  x/H  relative  to  unity;  that  is  x/II  =  x  <<:  1  ,  an 
approximation  that  is  discussed  in  detail  in  References  15  to  17  and  19. 
Thus  Eq.  (12)  may  be  written  for  small  time  as 


t 

i 

}  ; 

1 1 
\ 


30 


at 


(o) 


=  a 


o2s/o) 

_ z 

3x2 


(21) 


whore  Q!  -  7  /T  and  whore  Eq.  (19)  is  used  and  only  first-order 

c  h 

terms  are  retained.  Similarly,  Eq.  (15),  to  the  first  order  is 


39 


(o) 


ox 


-  t(7)fl)  +  i 


d‘) 


d  t 


(o) 


2  1 
V  "  V  +  ~l, 


(22) 
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and  the  first-order  parts  of  Eq.  (13)  and  (14)  are 


Io)(ji,0)  =  0 

(23) 

(r,)(0,b  .  b  <0>(t> 

V 

(24) 

The  reaction  time  is  neglected  since  it  is  generally  very  small 

until  S  becomes  substantially  greater  than  unity.  Now  write  Ea.  (22) 
v 

in  the  form 


where 


oo.  dS  (t) 

j6  A  V  A 

- -  (o  t)  -  A  -  =  -  B  t  , 

Ox  - 

dt 


(25) 


and 


3  =  (Y 


3/3  (rj  -  n  +  Y  l  x  ) 

(26) 

3(nti) 

(27) 

/  C  \  /  o  \  /  T  \ 

~  IT  - 

9 

(28) 

'  Cl'  '  °Z  '  '  *c  ' 

The  problem  consisting  of  the  differential  operation  (21)  with 

initial  conditions  (23)  and  boundary  conditions  (22)  and  (24)  belongs 

27 

to  a  class  of  problems  solved  by  Carslaw  and  Jaeger  using  Laplace 
transforms.  The  solution  for  small  time  is 
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V 


n  th 

Where  i  erfey  refers  to  the  n  integral  of  the  crfcy  function 
(see  Reference  27).  Consequently,  the  temperature  at  the  bubble 
surface  (x  =  0)  is 


„  (i) 


(0,  t) 


9(0)(t) 

v 


(30) 


In  comparison  with  the  temperature  in  the  bubble  for  an  adiabatic 
collapse  (no  heat  or  mass  transfer): 


(o) 


v 

a 

va 


7?  +  1 


T)  -  T?  + 


(31) 


whore  9  is  the  temperature  for  the  adiabatic  collapse  of  the  same 
bubble.  Generally  7)  is  in  the  range  of  10-20  so  the  initial 
temperature  rise  in  a  collapsing  vapor  bubble  with  heat  and  mass 
transfer  at  the  bubble  wall  may  be  almost  an  order  of  magnitude  less 
than  the  corresponding  rate  of  temperature  rise  for  the  adiabatic 
collapse  of  a  similar  bubble. 

The  temperature  profile  for  small  time  is  plotted  in  figure  1  for 
a  representative  set  of  parameters. 
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FIGURE  D-1  TEMPERATURE  PROFILES  IN  LIQUID  FOR  SMALL-TIME  SOLUTION  (EQ.  29) 
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ABSTRACT 


A  pair  of  piezoresistant  ytterbium  stress  gages  were  used  to 
record  pressure-time  histories  in  liquid  monopropellants  in  a  study 
of  the  shock  initiation  of  low  velocity  detonation  (LVD).  The 
propellant  was  contained  in  a  cylindrical  steel  tube,  and  the 
initiating  shock  was  produced  by  a  totryl  pellet-attenuator  donor 
system.  Variation  of  the  propellant,  the  wall  thickness-to- 
diameter  ratio  of  the  tubes,  and  the  initial  shock  pressure  from 
30  to  50  kbars  resulted  in  (a)  shock  decay  with  no  appreciable 
reaction,-  (b)  initiation  of  LVD,  and  (c)  initiation  of  high  velocity 
detonation  (HVD).  Peak  pressures  recorded  at  4  to  10  inches  down  the 
tube  were  in  the  1-  to  2-kbar  range  for  shock  decay,  and  in  the 
4-  to  7-kbar  range  for  LV1),  but  no  pressures  were  recorded  for  HVD. 
The  average  LVD  propagation  velocity  measured  with  the  gages  was 
about  2.0  mm/gscc  for  ethyl  nitrate  and  about  3.0  mm/gscc  for  FKFO. 
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Introduction 


«-  •"S'V  » 


This  paper  presents  an  experimental  investigation  of  the  shock 

initiation  of  low  velocity  detonation  (LVD)  in  liquid  monopropellants 

contained  in  cylindrical  steel  tubes.  The  object  of  the  work  was  to 

achieve  a  more  quantitative  description  of  LVD,  Since  LVD  propagation 

1,2 

velocities  had  been  measured  in  previous  studies,  and  the  circum- 

3 

ferential  strain  histories  had  also  been  examined,  a  better  description 
of  states  attained  in  the  liquid  was  required  to  achieve  this  objective. 
Piezoresistant  in-material  ytterbium  stress  gages  were  therefore  used 
in  the  experiments  to  measure  the  pressure  time  histories  of  liquid 
monopropellant  during  the  initiation  and  propagation  of  LVD.  A  pair 
of  gages  were  used  in  each  experiment  so  that  the  pressure  records 


Experimental  Procedure 


Figure  1  is  a  schematic  diagram  of  the  LVD  experiments  that  were 

performed  in  the  Stanford  Research  Institute  explosives  vault  facility. 

A  cylindrical  seamless  cold-rolled  steel  tube  was  used  to  contain  the 

liquid.  To  the  bottom  of  the  tube  was  attached  a  3-inch  diameter  disk 
* 

of  Homalite,  two  2-inch  diameter  by  1  inch  thick  pellets  of  tetryl, 
and  one  1/4-inch  thick  by  1/2-inch  diameter  tetryl  pellet.  The  tetryl 
detonation  was  initiated  by  an  exploding  bridgeware  detonator  in 
contact  with  the  small  pellet.  The  Homalite  disk  served  as  a  stress 
pulse  attenuator  and  its  thickness  was  chosen  to  produce  the  desired 
peak  pressure  in  the  liquid.  To  determine  this  pressure,  a  standard 
four- terminal  pie>  '.resistant  manganin  stress  gage  was  placed  in  the 
Homalite  1/32  o'  an  inch  from  the  surface  in  several  of  the  experiments 
to  record  ihe  peak  stress  that  propagates  into  the  liquid. 

Two  IT- shaped  4-terminal  piezoresistant  ytterbium  stress  gages, 
photoetchcd  from  0.001-inch  thick  foil,  were  mounted  on  a  long  strip 
of  0.003-inch-weave  fiberglass  cloth  whose  width  was  the  same  as  the 
inside  diameter  of  the  tube.  Long  thin  copper  leads  were  soldered  onto 
the  gage  terminals  and  a  layer  of  Homalite  was  spread  over  the  fiber¬ 
glass  to  contain  the  gages  and  leads  in  a  stiff  0.007-inch  thick 
package.  Fiberglass  fins  to  ensure  that  the  gage  package  remained  in 
the  center  of  the "tube  wore  attached  and  the  package  was  slid  into  the 
;ube  and  glued  into  place  with  the  leans  emerging  from  the  top. 

This  arrangement  of  the  gage  package  was  motivated  by  several 
factors.  First  of  all,  since  the  package  had  a  very  small  cross 
section  in  the  plane  of  the  expected  shock  front  (perpendicular  to 
the  tube  axis),  it  wouldn't  move  much  as  the  shock  and/or  detonation 
wave  jorcod  the  liquid  to  move  up  the  tube.  Second,  any  motion  of  the 
*  Trade  name,  Homalite  Corporation. 
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gage  package  would  not  stretch  and  break  the  leads  as  quickly  as  if  the 
leads  had  been  attached  to  the  tube  wall.  And  third,  since  the  coupling 
of  the  liquid  and  tube  wall  shock  responses  is  very  critical  to  the  LVI) 
phenomenon,  it  was  desirable  to  avoid  interference  with  the  integrity 
of  the  wall,  as  would  have  occurred  had  the  leads  needed  to  emerge  from 
the  walls.  One  disadvantage  of  this  orientation  of  the  gage  was  an 
increase  in  t.he  response  time  of  the  gage.  Since  the  thickness  of  the 
active  element  in  the  direction  of  shock  propagation  was  about  0.020  inch, 
the  fastest  rise  that  could  be  seen  by  the  gage  was  about  250  nsec. 

For  several  experiments,  one  or  two  Micro-Measurements  Type  EA-06- 

5006R-120*  strain  gages  were  glued  to  the  outside  of  the  tube  wall 

adjacent  to  the  location  of  the  stress  gages  on  the  inside  of  the  tube. 

They  were  oriented  to  measure  hoop  or  circumferential  strain  and  thus 

4 

to  compare  the  strain  histories  in  these  and  ea. Her  experiments, 
and  possibly  to  determine  the  time  lag  between  LVD  initiation  at  the 
center  of  the  tube  and  the  circumferential  expansion  of  the  tube. 

Eight  experiments  wore  performed  os  part  of  the  LVI)  program.  The 
relevant  experimental  parameter  for  each  of  the  shots  Is  given  in  the 
uppor  half  of  Table  1.  For  the  first  5  shots,  ethyl  nitrate  was  used 
for  the  liquid  monopropellnnt  and  wall  thickness-to-diamoter  ratios  and 
attenuator  thicknesses  were  chosen  that  had  previously  keen  successful 
in  LVI)  initiation  in  similar  materials.'  ’  ^  When  it  became  clear  that 
it  was  difficult  to  create  a  steady-state  LVI)  wave  in  ethyl  nitrate 
(see  next  section),  it  was  decided  to  switch  to  FEFO  (fluorodinitro- 

5 

ethylformal),  a  liquid  known  to  bo  more  sensitive  to  LVI)  initiation. 

After  shot  No.  6  failed  to  cause  detonation  in  the  FEFO,  the  liquid 

was  analyzed  and  was  shown  to  contain  a  considerable  amount  of  impurities. 

A  now  batch  of  FEFO  was  prepared  for  shot  Nos.  7  and  8. 

*  Micro-Measurements,  Romulus,  Michigan. 
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All  of  the  shots  were  instrumented  except  for  shot  No.  2,  whose 

purpose  was  simply  to  determine  by  terminal  observation  if  LVI)  had 

taken  place.  For  all  of  the  instrumented  shots,  a  pulsed  constant- 

current  power  supply  was  used  in  conjunction  with  the  ytterbium  and 

manganin  stress  gages  and  the  signals  were  recorded  on,  oscillographs  in 

6 

the  manner  described  by  Keough.  The  piezoresj s tance  coefficient  of 
ytterbium  (the  factor  relating  change  in  gage  'esistance  to  stress) 

7 

used  in  that  data  analysis  was  that  determined  by  Ginsberg. 
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Fxperimental  Kosults 


The  results  of  the  1AD  experiments  are  tabulated  in  the  bottom 
half  of  Table  1.  Terminal  observations  of  the  remains  of  the  steel 
tube  were  made  to  determine  what  type  of  detonation,  if  any,  the  liquid 
underwent.  Figaro  2  contains  photographs  of  the  tube  remains  for  the 
three  different  outcomes  that  wore  observed,  as  well  as  a  photograph  of 
the  tube  prior  to  detonation.  For  the  ease  in  which  no  liquid  deto¬ 
nation  took  place,  the  tube  is  undamaged  except  for  the  first  inch  or 
two  adjacent  to  the  explosive  pellet.  For  the  case  in  which  the  liquid 
underwent  I  AD,  the  section  of  tube  below  the  liquid  level  broke  into 
several  pieces  which  wore  often  severely  bent  or  twisted.  Some  of 
these  pieces  did  not  break  off  from  the  undamaged  upper  section  of  the 
tube  (the  section  above  the  liquid  lino)  but  bent  back  beyond  the  un¬ 
damaged  section  in  the  shape  of  a  partially  peeled  banana  (with  the 
fruit  removed).  'Ibis  was  particularly  true  of  the  ethyl  nitrate  shots 
No.  3,  '1,  or  f>,  where  analysis  of  the  gage  records  indicated  a  nonsteady 
state  I  AD,  as  will  be  discussed  shortly.  For  the  case  in  which  the 
liquid  underwent  UVD,  the  tube  broke  into  many  long  very  thin  fragments, 
as  would  be  expected  due  to  the  very  high  pressures  in  the  tube  during 
an  1 1 VI) .  For  shot  No.  7  in  which  II VI)  occurred,  no  stress  records  were 
obtained  duo  primarily  to  high  propagation  velocity  of  the  IIV1)  wave 
(between  6  and  7  mm/|iscc),  which  caused  the  wave  to  pass  the  ytterbium 
gage  location  long  before  the  triggering  of  the  gage  power  supplies  and 
oscilloscopes,  which  had  been  sot  to  expect  a  slowei  ••elocity  I  AD  wave. 

The  peak  pressures  recorded  by  the  various  stress  gage.'  are  also 
given  in  Table  1  and  the  average  propagation  velocity  of  the  peak  stress 
seen  by  the  ytterbium  gages  Tor  the  shots  In  which  both  gages  produced 
analy/.ablc  records.  A  typical  manganin  gage  record,  that  from  shot  No.  1, 
is  shown  in  Figure  3.  There  is  an  initial  peak  of  -v  f>9  kbar  which  is 


the  stress  that  first  reaches  the  gage  plane  in  the  Homalite,  followed 
by  a  drop  to  46  kbar,  which  is  the  stress  induced  into  the  liquid, 
and  is  lower  due  to  the  impedance  mismatch  between  the  Homalite  and  the 
liquid.  Since  the  attenuator  thickness  was  nearly  the  same  for  shots 
Nos.  1-5,  and  since  the  explosive  pellet  configuration  was  identical 
for  all  of  the  shots,  it  can  be  assumed  that  very  nearly  the  same 
stress  was  induced  into  the  ethyl  nitrate  for  shots  Nos.  1-5.  Manganin 
gage  records  for  the  shots  using  FKI'O  yielded  a  pressure  of  «s  43  kbar 
for  shots  Nos.  6  and  7  using  approximately  a  1-inch  thick  attenuator, 
and  ••  33  kbar  for  shot  No.  8,  which  used  a  lj-inch  thick  attenuator. 

The  ytterbium  gages  generally  worked  very  well  in  measuring  the 
liquid  pressure  histories  for  the  shots  in  which  the  liquid  underwent 
I  AD.  A  couple  of  the  ytterbium  records  (those  from  shot  No.  8) 
exhibited  a  significant  amount  of  high-froquonev  noise,  but  the 
remainder  of  the  records  had  excellent  signal- to-noisc  ratios.  The 
records  from  the  two  ytterbium  gages  in  shots  Nos.  3  and  5  are  shown 
in  Figures  4  and  5,  respectively.  Observation  of  the  various  records 
from  shots  Nos.  3,  4,  and  5,  however,  showed  that  the  pressure  histories 
at  the  two  locations  are  qualitatively  dissimilar  for  each  shot.  This 
is  most  clearly  seen  in  Figure  4.  The  first  gage,  located  at  a  distance 
of  6  inches  from  the  detonating  end  of  the  tube,  recorded  a  rapid  rise 
(less  than  1  psoc)  to  a  peak  pressure  of  about  2^  kbar  and  then  a  slow 
decay  to  zero  pressure.  Furthermore,  the  gage  remained  continuous  for 
at  least  20  nsec  following  the  initial  rise.  On  the  other  hand,  the 
second  gage,  located  at  10  inches  from  the  end  of  the  Lube,  exhibited 
a  slower  rise  (more  than  3  usee)  to  a  stress  of  over  5  kbar  and  then  a 
nearly  level  plateau  until  the  gage  began  to  stretch  and  break  3  nsec 
later.  The  latter  gage  record  is  more  like  that  to  be  expected  in  an 
1AD,  while  the  lower  stress  and  longevity  of  the  former  is  indicative 
of  no  I, VI)  in  the  immediate  vicinity  of  the  gage  package.  The  records 


from  shot  .Vo.  5  exhibit  similar  tendencies,  as  shown  in  Figure  5, 
except  that  here  the  gage,  located  'lA  inches  from  the  end  of  the  tube, 
has  a  higher  peak  pressure  (4.1  kbar)  and  shorter  recording  life  than 
the  gage  located  7A  inches  from  the  ends  (with  a  peak  pressure  of  2.7 
kbar).  Since  previous  experimental  studies  and  computer  simulation 
indicate  that  I. VI)  should  start  at  distances  c-f  less  than  a  few  inside 
diameters  down  the  tube,  it  is  reasonable  to  assume  that  if  a  steady- 
state  I.VI)  were  to  bo  initiated  in  the  tube,  it  would  begin  before  the 
pressure  pulse  reaches  the  first  gage  location.  Therefore  it  appears 
that  a  steady-state  I  AT)  wave  was  not  produced  in  the  region  be  tween  the 
two  ytterbium  gages.  More  likely  a  pulsating  or  oscillatory  I, VI)  was 
produced,  or  perhaps  merely  a  few  isolated  hot  spots  wore  initiated  at 
various  points  in  the  tube.  The  higher  pressure  seen  by  one  of  the 
gages  in  both  shot  No.  3  and  No.  5  as  compared  to  shot  Vo.  4,  the 
higher  peak  pressure  propagation  velocity  (about  2.00  mm/psoc  as  compared 
to  1.50  nun/ gso c) ,  and  the  fact  that  the  tube  in  shot  No.  4  was  less 
damaged  than  the  tube  in  shots  3  and  5  indicate  that  significantly  less 
of  the  liquid  in  shot  No.  4  underwent  I.VI)  than  in  shots  No.  3  or  No.  5. 

Many  explanations  are  possible  for  the  failure  to  Initiate  a  steady- 
state  l.VI)  in  ethyl  nitrate.  Insufficient  experimental  data  are 
available  to  predict  precisely  the  conditions  under  which  such  .a  phenom¬ 
enon  will  occur.  It  is  therefore  very  possible  that  the  experimental 
parameters  chosen,  such  as  the  wall  thicknoss-to-tube  diameter  ratio 
and  attenuator  thickness,  were  slightly  different  from  those  needed  to 
achieve  steady-state  1AD.  Perhaps  ethyl  nitrate  is  sufficiently  insensi¬ 
tive  in  the  region  of  those  parameters. 

The  two  stress  gage  histories  for  shot  No.  8  using  FHFO  as  the 
liquid  were  much  more  qualitatively  alike,  although  high  frequency 
noise  obscured  their  structure  somewhat.  Furthermore  the  higher 
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measured  peak  pressures  (4  to  7  kbar)  and  faster  average  propagation 
velocity  (3  mm/psec),  despite  the  liquid  being  hit  by  a  lower  stress, 
indicated  that  a  steady-state  LVI)  was  more  likely  initiated  in  FEFO 
than  in  the  ethyl  nitrate. 

The  final  experimental  result  to  be  reported  was  that  of  the 
strain  recorded  for  shot  No.  4  (the  strain  gage  records  in  shot  No.  1 
were  too  noisy  to  analyze).  The  peak  circumferential  strain  seen  at 
the  distance  of  10  inches  from  the  end  of  the  tube  was  0.2%.  Comparison 
with  the  computer  simulation  reported  in  Appendix  A  indicated  that  a 
0.2%  strain  corresponds  to  an  internal  liquid  pressure  of  a  few  kilobars, 
which  is  just  that  recorded  by  the  stress  gage.  Furthermore  LVI)  would 
produce  a  circumferential  strain  in  the  tube  wall  of  more  than  1% 
shortly  after  the  wave  passed  that  section  of  the  tube.  This  is  yet 
another  indication  of  the  absence  of  LVI)  in  the  region  10  inches  from 
the  detonated  end  of  the  tube  for  shot  No.  4. 
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Conclusions 

4 

Although  peak  pressure  measurements  have  been  reported  the  results 

reported  here  are  the  first  direct  measurements  of  pressure  history  in 

a  liquid  undergoing  LVI).  In  the  shots  in  which  successful  LVI)  was 

initiated,  peak  pressures  of  about  4  to  7  kbar  were  observed  which  is 

1-4 

consistent  with  other  observations  of  LVI).  The  measurements  were 
taken  with  a  minimum  of  interference  with  the  detonation  wave  by  the 
use  of  very  small  gages.  Therefore,  the  rise  times  and  pressure 
history  measured  by  the  gages  are  believed  to  be  closely  representative 
of  the  true  stress  histories  in  an  LVI)  wave.  This  technique  should 
prove  useful  for  recording  pressure  history  ahead  of  and  within  LVI) 
waves  for  comparison  with  theoretical  prediction  and  to  gain  a  better 
understanding  of  LVI)  behavior  in  various  liquids  as  a  function  of 
confinement,  initiating  shock  strength,  and  other  experimental  parameters. 
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FIGURE  E-3  OSCILLOGRAM  FROM  MANGANIN  STRESS  GAGE 
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FIGURE  C-l  OSCILLOGRAMS  FROM  SHOT  HO  3  FOR  THE 
Vn FRRIUM  STRESS  GAGES  AT  DISTANCES 
SHOWN  IN  (,  1  AND  lb) 
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FIGURE  E-5  OSCILLOGRAMS  FROM  SHOT  NO.  5  FOR  THE 
YTTERBIUM  STRESS  GAGES  AT  DISTANCES 
SHOWN  IN  (a)  AND  (l>) 
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